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At the request of the Eev. Charles Cheyne, the father of the 
late C. H. H. Cheyne the author of this work, who to the 
great regret of all his relatives and friends died at Torquay 
on the 1st of January, 1877, 1 have undertaken the charge 
of this new edition. I have made but two changes in the 
body of the work, and have introduced but few notes. Inas- 
much as the methods of Chapter II. are peculiarly the 
author s own I have not ventured to alter them : but I have 
instead added to the Appendix the Corollary at the end of 
Art. 5, a method of determining the variation of the elements 
(Art. 6) which is for the most part due to Pont^coulant with 
a slight addition of my own to make it complete, together 
with Art. 8 also largely due to the same author. I have 
also given some additional examples and have indicated the 
sources from which they were obtained. 

If a new treatise on the Planetary Theory had to be 
written it would be necessary to digest the work of Jacobi, 
Hansen, Leverrier and Newcombe, to mention only a few 
eminent analysts. A brief account of a variety of methods 
might be found in Matthieu's Dynamique. 
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VI PREFACE. 

A biographical notice of the late author may be read in 
the Monthly Notices of the Royal Astronomical Society, 
Vol. XXXVII. pp. 147, 148. This I find as complete as is 
required in such a sketch of the too brief life of a diligent 
lover of Mathematical and Astronomical Scienca 

A. FREEMAN. 



MuBSTON Bectory, Kent, 
October 31«t, 1883. 
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PREFACE TO THE FIRST EDITION. 



In this volume, an attempt has been made to produce a 
Treatise on the Planetary Theory, which, being elementary 
in character, should be so far complete, as to contain all that 
is usually required by students in this University. But it 
is not without diflSdence that I submit my volume to their, 
noticed In the earlier part of it, the methods which have 
been adopted are to some extent original*, and the general 
arrangement of the second Chapter will, it is believed, be 
found to be new. Through the kindness of the Publishers, 
a portion of Pratt's Mechanical Philosophy has been placed 
at my disposal. Of this I have availed myself, particularly 
in the Chapter on the Stability of the Planetary System; but, 
on the whole, comparatively little has been reprinted ver- 
batim from that work. Among other sources of information, 
my obligations are mainly due to Pontdcoulant's Th^orie 
Analytique du Systime du Monde, Airy's Mathematical 
Tracts, and Frost's Planetary Theory in the Quarterly 
Journal of Mathematics : but I have also referred to the 

* Some of these have already appeared in Mathematical Journals. 
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Micanique Cileste^ the Micanique Ancdytique^ Mrs Somer- 
ville's Mechanism of the Heavens, (a work forming a complete 
Mathematical Treatise on Physical Astronomy^) a Memoir 
by Prof. Donkin on the Differential JSqiuUions of Dynamics, 
Phil. Trans, 1855, &c. A collection of Problems has been 
added, taken chiefly from the Smith's Prize and Senate- 
House Examination Papers of the last twenty years. In 
conclusion, 1 would express my sincere thanks to Messrs. 
A. Freeman, P; T. Main, and other friends, of St John's 
College, for the valuable assistance which they have afforded 
me, and would venture to hope that the work will be found 

useful. 

C. H. H. CHEYNE. 

St John's College, 

October, 1862. 



In the Second Edition comparatively few changes have 
been made. The work has been revised, and, it is hoped, 
in some degree improved. The Stability of the Planetary 
System has been rather more fully treated, and an elegant 
geometrical explanation of the formulae for the secular 
variations of the node and inclination introduced, for which 
I am indebted to a paper by Mr H. M. Taylor, Fellow of 
Trinity College, in the Oxford, Cambridge and Dublin 
Messenger of Mathematics, 

C. H. H. C. 

1, Dean's Yabd, Wesxminsteb, 
September, 1870. 
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CHAPTER I. 



INTRODUCTION. 



1. To determine the motion of a system of bodies sub- 
ject only to their mutual attractions, is a problem the mathe* 
matical difficulties of which have not yet been overcome: 
hence in the particular cases of this problem which Physical 
Astronomy presents, recourse must be had to methods of ap- 
proximation. Happily the arrangement of the Solar System 
renders approximate methods possible, and in the skilful 
hands of the Mathematicians of the last century, they have 
been brought to a high state of perfection. 

2. If the Sun were the only attracting body, the planets 
would describe exact ellipses, agreeably to Kepler's first law ; 
but in consequence of the attractions of the planets them- 
selves, slight deviations from elliptic motion are produced. 
The method of calculating these deviations, to which our 
attention will chiefly be directed, is due to Euler, but was 
subsequently extended and perfected by Lagrange : it con- 
sists in supposing the planets to move in ellipses, the ele- 
ments (or arbitrary constants) of which are continually though 
slowly changing*. 

* The legitimacy of this hypothesis wiU appear when we come to treat 
of the equations of motion. See Arts. 21 and 22. 

C. P. T. 1 
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2 PLANETARY THEORY. 

3. Now the elements of an elliptic orbit are (i) the mean 
distance, or semi-axis major, (ii) the excentridty, (iii) the 
longitude of perihelion, i.e. of the point of the orbit nearest 
to the Sun, (iv) the longitude of the epoch*, or mean longi- 
tude at the epoch from which the time is reckoned, (v) the 
inclination of the plane in which the orbit lies to some fixed 
plane of reference, (vi) the longitude of the ascending node. 
Of these (i) and (ii) determine the magnitude of the orbit, 
(iii) determines its position in its own plane, (v) and (vi) de- 
termine the position of this plane, and (iv) has reference to 
the position of the body itself in its orbit. 

If the planets moved accurately in ellipses, these would be 
constants : we must however be prepared to consider them as 
variable quantities, which it will be the object of the problem 
to determine. They are termed the elements of the orbit 

4. But further, npt only is it found that the true orbit 
of a planet is not an ellipse, but that it is not even a plane 
curve, although the departure of the planet from the plane 
in which it is at any instant moving is extremely slow. We 
define as the plane of the orbit the plane containing the 
radius vector and direction of motion of the planet at the 
instant under consideration. 

5. We shall suppose the Sun and planets so distant from 
each other that they may be considered to attract as if they 
were condensed into their respective centres of gravity; a sup- 
position which would be rigorously true if these bodies were 
exactly spherical, and either of uniform density or composed of 
concentric spherical shells, the density of each shell being uni- 
form throughout. The errors, however, thus introduced into 
the motions of translation are found to be inappreciable for the 
planets, though not in the case of their satellites. The mo- 
tions of rotation will not be considered in the present treatise. 

* Also briefly termed the epoch. 
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6. Moreover, since the masses of the planets are ex- 
tremely small in comparison of that of the Sun, it follows 
that in cases where it is not necessary to carry the approxi- 
mation beyond the first order of these masses, we are permit- 
ted to avail ourselves of the Principle of the Superposition of 
Small Motions, and thus to reduce the problem to a case of 
that of the Three Bodies. 

7. So far the Theory of the Planets resembles that of 
the Moon, and the same method of treatment might be em- 
ployed in both cases. But they differ in this respect: the 
ratio of the distances of the disturbed and disturbing bodies 
from the central one* is much smaller in the Lunar than in 
the Planetary Theory, so that if in the latter theory the ap- 
proximation were made by means of series proceeding by 
powers of this ratio, it would be necessary to retain many 
more terms than are required in the former. On the other 
hand, the perturbations of the Moon are far larger than those 
of the planets, since in the former case the Sun, of which 
the mass is enormous, and the distance not proportionately 
great, is one of the disturbing bodies. For these reasons 
different methods of calculation are employed. 

8. To find an expression for the component in any di- 
rection of the force which disturbs the motion of a given planet 
relatively to the Sun, 

Let M denote the mass of the Sun, m, m', m"', &c., those 
of the planets, and suppose the relative motion of m re- 
quired. 

Let Xy y, z, x, y\ /, a?", y", /', &c., be the co-ordinates 
of w, m', m\ &c., referred to any system of rectangular axes 

* By the central body is meant that whose attraction exercises the 
greatest influence on the body whose motion is required; the Sun, for 
instance, in the Theory of the Planets, and the Earth in that of the Moon., 
All the other attracting bodies are called disturbing bodies. 

1—2 
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originating in the centre of gravity of the Sun ; r, /, r", &c., 
their distances from the origin ; f! ^ p\ &c., the distances of 
m^ m'^j &c., from m. 




Now if to every body of the system we apply forces equal 
and opposite to those which act upon the Sun, we shall re- 
duce the latter to rest without aflfecting the relative motion. 
Hence, considering the action of only one disturbbg planet 
m', the forces acting upon m will be 

— -^ — in direction mM, or briefly, ^y 

-7« in direction mw!. 
P 

m' . 

-75 in direction m'M, 

r 

of which the last two constitute the disturbing force. 

On the hypothesis of Art. 5, the components parallel to 
the axes of the forces acting on m, will be 



r p r 

^fMZ mf (z — z) mfz 

of which the terms containing rn are the disturbing forces. 
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Let 8 denote the length of the arc of any curve measured 
from some fixed point up to m: then the resolved part of the 
disturbing force parallel to the tangent at m to this curve 
will be 

mf {x' — x) dx ifn! {j/ ^ y) dy w! (/ — 5?) dz 

^ fmx' dx m'y' dy mz' dz\ 
^\l^d8^'7^"d^'^~^ dsj' 



d (m' iri , , , -.) 



ds 

on supposition that Xy y, Zy are alone aflfected by the process 

dR 
of differentiation, which may be written —^ , if 

If we express in like manner the disturbing forces due 
to the action of m", w"", &c., we shall have for the whole 
component in this direction 

d^ d_^ 
ds da '"' 

or^,ifi2 = i2'+i2" + .... 
ds 

The function It is called the disturbing fu/nction: the 
arbitrary curve employed above will be termed the curve of 
reference. 

9. From the manner in which -y- has been introduced, 

ds 

it appears that R is supposed to be expressed in terms of s 
and quantities which do not vary with s. It must, however, 

be borne in mind that in -r- the variation is purely hypo- 
thetical, and has nothing whatever to do with the actual 
variation of jB due to the motion of the planet. 
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For example, suppose the curve of reference a straight 
line parallel to the axis of Xy and let R be expressed in 
terms of a?, y and z ; then in this case x only will vary, and 
the disturbing force pai*allel to the axis of x will be denoted 

by -T- , y and z being considered constant in the diflferen- 

tiation. Similarly, the disturbing forces parallel to the axes 

of y and z will be expressed by -r- and -7- respectively, the 

diflferential coefficients being strictly partial. 

Again, suppose the cur\re of reference a circle with its 
plane parallel to that of xy, and its centre in the axis of Zy 
and let R be expressed in terms of the polar co-ordinates 
(rj, 6^ of the projection of the planet on the plane of xy, and 
its distance {z) from this plane ; then in this case 0^ only will 
vary, and the disturbing force perpendicular to the projected 

radius vector will be expressed by —-13- , r^ and z being con- 
sidered constant in the diflferentiation. Similarly, the forces 
parallel to the projected radius vector and to the axis of z, 

will be expressed by the partial differential coefficients -j- , 

dR .. 1 

-T- respectively. 

10, The disturbing function is independent of any par- 
ticular system of co-ordinates that may be employed. For 

r 



p r^ \r r' r r* rr) 



9 / 
m mr 

-7- - --TT cos ft), 
p T 



if to denote the. inclination of r to r. 
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11. To express R' in terms of the polar co-ordinates of 
the projections of m and m' on a fixed plane, and of their 
distances from it 

Take the fixed plane for that of xy : let r^, r/ be the pro- 
jections of r, r' upon it, and 0^, 6^ the inclinations of r^, r^ to 
the axis of x ; then 

a;=*rjCOS^j, y=^ r^AnO^, 

X = r/ cos ^/, y' = r/ sin ^/; 

therefore awe' + yy' + W = r// cos (^^ — ^/) + zz\ 

p'^^[x^xy+{y^yy^{z^zy 

= r,' + r;'-2r,r;cos(^,-^/) + (^-/)', 

r' =^'* + y'» + ^;'« 

Hence by substitution, 



h' + »•.'* - 2Vi' cos (^, - 0{) + {z- «•)'}* 

ot' {r.r/ cos (0^ - 0^') + zz\ 

12. In a subsequent Chapter we shall consider the de- 
velopment of U in terms of the time and the elements of the 
orbit, in a series ascending by powers and products of the 
eccentricities and inclinations, which for the principal planets 
are very small*. At present we shall content ourselves with 
shewing h)w R may be expressed in terms of these quantities. 

* For the smaller planets and comets this is not the case; so that 
different methods of calculation are required for these bodies. 
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We shall assume that the equations connecting the co-ordi- 
nates^ the time, and the elements in an elliptic orbit, hold 
also when the motion is disturbed. 



13. To explain how R may be expressed in terms of Hie 
time and the elements of the orbit 

Let r, 6 denote the radius vector and longitude of the 
disturbed planet, the latter being measured on a fixed plane 
of reference as far as the node, and thence on the plane of the 
orbit : let the elements be a the mean distance, e the excen-. 
tricity, w the longitude of perihelion, e the longitude of the 
epoch, (the last two being measured in the same way as ^,) 
A the longitude of the node measured on the plane of refer- 
ence, and i the inclination of the plane of the orbit to the 
plane of reference. Our object is to express R in terms of t 
and these elements. 

Again, let 0^, m^, €^, Sl^ denote the longitudes of the 
planet, of perihelion, of the epoch, and of the node, measured 
entirely on the plane of the orbit. 

Let a sphere be described with its centre coinciding with 
that pf the Sun, and its radius of any magnitude : let the 
planes of reference and of the orbit cut it in the great circles 
NM, NPy then the lines of nodes will cut it in N] let the 
radius vector of the planet cut it in P, the projection of this 
radius on the plane of reference in Jf, and the lines from 
which By Oq are measured in i, respectively. We shall 
suppose L to be the same origin as that from which ^^ is 
measured in Art. 11. 

Then in the figure LM=0^, LN' + NP = 0, OP = 0,, 
LN=ily the angle FNM=i, and Pif = the latitude of the 
planet which we shall denote by X. 
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Hence from the right-angled triangle PNM, 

tan(5j -n) = cost tan (^-fl) (1), 

sin\ = sint sin(^ — fl) (2): 

also rjS=r cosX (3), 

z=^r sinX (4). 




Again, from the formulaB of elliptic motion*, 

r = a {1 + Je' — e cos (n< + €<, — -cr J — Je* cos 2 (n^ + e^ — tsr^) — . . . } , 

0^ = nt + €^ + 2eBm{nt + €^-'ar^)-\'^ sm2{nt + €^'''GTo)+..,f: 

but ^-^, = iJ\r- 0J\r = €-6, = t!r-t;T,; 

therefore O^O^ + e — e^, €f^-'vr^ = € — w, 

and our formulae become 

r=a{H-Je'—CCOs(7i<+€—'sr)—J^e*co82(nit+e— tar)— ...)... (5), 

0=int + € + 2e 8m{nt+ €- Iff) +|e'sin2(n^+e-t!r)+ (6). 

In Art. 11 we have expressed jR' in terms of r,, 0^ and z ; 
hence by equations (1) to (4) it may be expressed as a func- 
tion of r, 0, fl, and i : w6 may then substitute for r and 
from equations (5) and (6), and jB' will be expressed in terms 

* See Tait and Steele's Dynamics, Art. 162, Fifth Editi<m. 
f n is termed the mean motion, and is conneoted with the mean distance 
by the equation n'a'=/*, where /*=ilf+»iik 
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of t and the elements of the orbit. Similarly iJ", JB'", &c., 
and therefore It may be expressed in terms of t and the 
elements. 

14. We proceed to investigate certain relations which 
subsist between the partial differential coefficients of B with 
respect to the co-ordinates of the disturbed planet, and its 
partial differential coefficients with respect to the elements 
of the orbit. These will be useful in obtaining the formulae 
by which the values of the elements are calculated. 

We premise that when we speak of the partial differ- 
ential coefficient of jB with respect to one of the elements, 
we suppose R expressed in the manner indicated in the last 
Article, and that the time as well as the other elements are 
considered consts^nt in the differentiation : when we speak of 
the partial differential coefficient of jR with respect to r qt 0, 
we suppose R expressed in terms of r, 0, i and il, which may 
be done by equations (1) to (4). 

15. To shew that t^ = -3- + -t— . 

d^ ae dts* 

Equations (5) and (6) of Art. 13 may be written 

^ — w = ^ (w^ + € — w), 
whence it follows that 

de d'sr ' 

d0 d0 __ - 
de dm 

Now since c and «r enter into JB only through r and 0y 

dR^dR dr dR d0 
de " dr de d0 de ' 
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dR^dR dr dR d0 ^ 
d'sr dr d«r dd dw ' 

therefore, by addition, 

dR dR ___ dR 
de dm dO * 

IP rp . ^. ^dR dR . dR dR 
lo. lo shew that ttt = -j-^ + j — f- tfs • 

d^j de dtiT dil 

From equations (1) to (4) of Art. 13 we obtain 

where 0, ;^, ^^ are symbols of functionality. 
It follows that 

^i + ^-O 
dd ^ dil "■ ' 

t/^ ''dil" ' 
dz dz ^ ^ 

dd'^da" 

Now since by Art. 11, i? is a function of r^, 0^, and z, 

dR^dR dr\ dR dO^ dR dz^ 
dd dr^ de^ dd^ dd^ dz dd' 

dR ^dR dr dR d0 dR dz^ ^ 
dil dr^ dn "^ d0. dil '^ dz dCl' 



therefore, by addition, 



dR dR __ dR 
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whence, by the last article, 

dR ___ dR dR dR 
dl^j "" de dvT dH ' 

17. To obtain , - in terms of -^ and ,^ . 

If w denote the excentric anomaly, we have* 

r = a(l — 6 cosw) (1), 

'«"4^-yG^)-i (^'' 

n< + € — -BTsi^ — ^ sinw (3), 

from which r and ^ may be expressed in terms of t and the 
elements by eliminating u. Assuming r and 6 so expressed, 

we proceed to obtain -y- and -r- . 

de de 

dr f . du 



From (1), J- = ct f e sinw -^ — cosw j , 



du 
and from (3), -y- (1 — ^ cosw) — sinw = (4) ; 

diL 
eliminating -j- , we have 

dr ( e sin*u 

jT- = a •{ q cosM 

flte [1 — e cosu 



_ e — cos-w I 
"" l — e cosu) ' 



(1 + e) sin' 2 - (1 - e) cos" ^ 



(1 - e) cos' 2 + (1 + ^) sin' ^ 

* See Tait and Steele's Dynamics^ Arts. 160 and 161, Fifth Edition, 
This nse of the auziliaiy u was first suggested by Mr W. Pirie, Carnb, Math. 
Journal^ 1837, Vol. i., p. 47. 
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tan«^-l 



^ a 



1 + tan' 



,<?~tBr 



, by (2), 



= — a cos (^ — -cj). 
Again, dififerentiatlng the logarithms of equation (2), 

2 \1 + e "^ 1 - J "^ si] 



sin {0 — zct) de 
eliminating -=- by means of (4)3 

1 d0 1 



sin u de ^ 



+ 



Bm{0'-"Gr) de l—e* l — ecost^' 

if h^ = /Aa (1 — e*) ; therefore 

Now since JR is a function of e only because it is a function 
of r and 0, 

dR^dRdr dRd0 
de ~ dr de d0 de ' 

Since 5 — -07=5^ — -07^, (see Art. 13,) this equation may 
be written 

under which form it will be useful in the next Chapter. 
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FOBMUL^ FOR CALCULATING THE ELEMENTS OF THE ORBIT. 

18. We now proceed to form equations of motion, 
taking the Sun's centre for the origin of co-ordinates, the 
radius vector of the planet for the axis of ocy a perpendicular 
to it in the plane of the orbit for the axis of y, and a normal 
to this plane for the axis of z. With this system, it will 
be shewn that two of the resulting equations can be ex- 
pressed in the same forms as if the planet moved in one 
plane. 

Let X, y, z be the co-ordinates, w, v, w the velocities of 
the planet with reference to three rectangular axes originat- 
ing in the Sun's centre, and moving with angular velocities 
^i» ^2» ^3 9'bout their instantaneous positions : let X, Y", ^be 
the accelerations due to the impressed forces in the directions 
of the axes. Then (Routh's Rigid Dynamics^ Arts. 244, 245 
Third Edition)^ 

doc , . 



V = -^.- — z(f)^+ or^j 



dt 



^^■£-^<l>2+y4>t 



(1), 



and the equations of motion are 
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,(2). 



^ dw ... 

In these equations ^j, (^^, (f>^ are arbitrary; we propose 
so to determine them that the axis of x may coincide with 
the radius vector of the planet, and the plane of xy with the 
plane of the orbit. 

In order that the axis of w may coincide with the radius 
vector of the planet, we must have 

a? = r, y = 0, js ^0, 

always ; and therefore 

dx^dr dy_ dz_ 
dt^dt' dt ' dt^ ' 

and in order that the plane of xy may coincide with the plane 
of the orbit, we must have 

always ; and therefore 

dt "• 

Hence equations (1) give 

dv 

^=di' *=''^»' i>^=^> 

and equations (2) become 



dt 



dt' 
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Now let <^3 = -v^ : thus for the equations of motion 
we have 



dt) ' 



rdt\ dtj' 

of which the first two are the same in form as if the plane 
pf the orbit were at rest. 

19. If we measure on this plane from the planet*s 
radius vector in a direction contrary to that of motion 
an angle equal to 0^, we arrive at what may be considered 
as the origin from which 0^ is measured. Since this' will be 
a point having no angular velocity about the axis of z, which 
is normal to the plane of the orbit, it is said to be fixed in 
the plane of the orbit*, 

20. We shall for the present confine our attention to the 
first two of the above equations. 

In order to find the components of the disturbing force 
parallel and perpendicular to the radius vector of the planet, 
let us take first the radius vector as the curve of reference ; 
then s = r, and R being supposed expressed as a function of 
r, 0, il, and i (see Art. 13), we have 

dR^dRdr dRd0 dRdn dR di 
ds dr ds d0 ds dil ds di ds 

_dR 
^ dr' 

* It is necessary to define the expression /aj^d in the plane of the orbit, 
since the definition of Art. 4 is not sufficient completely to regnlate the 
motion of this plane, though affording it a distinct geometrical position. 
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since 6y XI and i do not vary with «, Henee the disturbing 

dR 



force in direction of the radius vector = 



dr " 



Again, let us take as the curve of reference a circle in the 
plane of the orbit, with its centre coinciding with that of the 
Sun ; then Bs = rS0, and we have 

dS^ldR 

ds ~r dd' 

since r, ft, and i do not vary with s. Hence the disturbing 
force perpendicular to the radius vector = - -^ . 

We have then 

y, __ fjb dR v" — ^ ^^ 

and the equations become 

dt\ \dt) 1*^ dr ^^'' 

^ /_. ddy dR .„ 

dtK~di)-'dd ^^^' 

21. These equations do not admit of rigorous integra- 
tion, but we may reduce them by the method of the Varia- 
tion of Parameters to a system of differential equations of the 
first order. The principle of this method may be explained 
as follows. Suppose it required to integrate the equations 



4>i \^y Vy ^y 



.(i), 



dx dy d^x d?y\ _ p 

~dt' di' de' dej" ' 

*P»r^'^'dr dt' df df) ^« 

where P^, P, are functions of t The solution of these equa- 
tions can be made to depend upon that of the equations 

C. P. T. 2 
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<f>^ = 0, ^g = 0. Suppose the four first integrals of these last 
equations to be 



(dx dv\ 

'''V'^'dt' dt) = '^ 

(dx dy\ 

'"•y'^'df dt)'='' 

( ^ dx dy\ 



(dx dy\ 
"'y'^'-dt' dtr'^ . 



(ii). 



where c^, c,, Cg, c^ are arbitrary constants or parameters. 
The method of the Variation of Parameters consists in so 
determining c^, c^, c, and c^ as functions of t, that these inte- 
grals (and therefore the two final integrals of the equations 

^^ = 0, <^g = 0, which can be obtained from equations (ii) by 

dx dij\ 

eliminating -j- and -^1 shall satisfy equations (i). That 

^1, c,, Cj, and c^ can be so determined, may be seen as follows: 

by tie solution of equations (i), values of x and y and there- 

dx dv 
fore of -t: and -rr can be found as functions of t and constant 
dt dt 

quantities ; if these be substituted in equations (ii) the requi- 
site values of c^, c^, Cg and c^ will be obtained. For an ex- 
ample of the application of this method, see Boole's Differen- 
tial Equations, Chap. IX. Art. 11. 

22. If in equations (1) and (2) of Art. 20 we put iJ = 0, 
and then integrate them, we obtain 



- = ^,{l+<?cos(^,-'cr,) 



dr fie . ,^ \ 



(3), 



(4). 
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^^-f^ (')' 

where h, e, tsr^ are the constants of integration. Equation (3) 
indicates motion in an ellipse, of which e is the excentricity, 
'sr^ the longitude of perihelion, and h twice the area described 
in an unit of time. If the mean distance in this ellipse be 
denoted by a, we have in addition 

A« = /Aa(l-6') (6). 

We shall assume (in accordance with the principles of the 
method of the Variation of Parameters) the first and second 
integrals of equations (1) and (2), together with equation (6), 
to retain the same forms when M is not zero ; h, e, ^q, and a 
being in this case considered variable*. 

The values of these elements are to be obtained from the 
condition that the above integrals shall satisfy equations (1) 
and (2). 

23. If their values as calculated for any given time be 

substituted in equation (3), it will represent an ellipse having 

a contact of the first order with the actual orbit, since the 

dr dO 
values of -r: and -rr^ at the common point will be the same 
at dt ^ 

for both curves. It is termed the instantaneous ellipse, since 

the planet may for an infinitely small time be supposed to 

move in it. Moreover, the velocity and direction of motion 

of the planet will be the same as if it moved in this ellipse, 

so that if at any time the disturbing force were to cease, the 

planet would continue to move in the instantaneous ellipse 

constructed for that time. This is accordingly sometimes 

given as the defijiition of the instantaneous ellipse. 

* We shall also for convenience suppose the equation n^a*=i/i to hold in 
the disturbed orbit, n being of course considered yariable. 

2—2 
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24. To ohtain formulce f(yr calculating the elements of the 
instantaneous ellipse at any time. 

Suppose the value of c required, where c denotes any one 

of the elements. From equations (3), (4»), and (6) we may, 

by eliminating the other elements, obtain o as a function of 

dv 
r, 6t,y -T7 , and A*: let then 
^ at 

c=f(r,0,,r\hl 

dv 
where i^ is written for tt • Differentiating, we have 

da^^dr df dd^ df_ d^ dfdh 
Jt~ drdt dd^ dt di^ dt dh dt' 



dr' ^r 



i'di) ~ S ■*■ d^ ' ^^"^ equation (I) ; 

dfdBdfdR 
^ dr' dr ^ dh dd' 

But, since by hypothesis, if B were zero and c constant, 
our assumed integrals would still satisfy the dififerential equa- 
tions, we have (making R zero and c constant) 

r,_dfdr d^de, df[(de,\'_!l^^ 

" ~ drdt ■*■ de, dt ^ dr' \\dt) r')' 

,^ , do df dR , dfdR 

therefore - = ^_+^_. 

* We retain h for convenience^ in preference to replacing it by r' ~ . 
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Hence in obtaining the formulaB for calculating the ele- 
ments of the orbit, we may proceed as follows. From equa- 
tions (3), (4), and (6) we may express the element required 

as a function o{ r, 0^, -j- , and h. We may then, by differen- 
tiating the resulting equation with respect to ^ as if r and 0^ 

, , ... dR « d^r J dR « dh , 

were constants, writing -p for -^ , and -jn ^^^ "ii * *^^ 

eliminating, if necessary. | and A by means of equations 

(4) and (6), obtain the differential coefficient of the element 
required in terms of the elements, the co-ordinates of the 
planet, and the disturbing force. The result, however, will 
in every case admit of being expressed in terms of the ele- 
ments and of the differential coefficients of R with respect to 
them, a form under which it is very convenient of applica- 
tion. 

26. To obtain a formula for calculating the mean dis^ 
tance. 

From equations (3), (4) and (6), if e and -cr^, be eliminated, 
we shall find 



a "~ \3i) r* r 



Differentiating as if r were constant, and writing -f- for 
cPr dR r dh , 

a'dt" [di^di'^dS'di)' 
Now since r =/ (nt + e — «r), 

^^ — 'Bjjj = ^ — -BT = ^ (w< + e — bt), 

tJv* /7/9 

and the forms of -r: and -rr are the same as if the elements 

dt dt 

were constant, we have 
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and similarly. _.= «_. = „_; 

a* dt ~ \dr de dB de/ 

da ^2na^dR 
dt fi de ' 

26. This formula may also be obtained as follows. If 8 
denote an arc of the actual path of the planet measured from 
some fixed point to its position at time t, we have the equa- 
tion of motion 

d^8 /* ^^ I ^ 

and by a known formula 

/d8\* ^ 2/A _ /A 
\dt/ r a * 

Differentiating the latter we obtain 

^d»d^ 2fidr jida 

Jtde" Vdt'^a^'dt' 

da 
and, multiplying the former by 2 -i- , 

^dsd^8_ 2fidr ^dRds 

^,1 u da c% dR ds 

therefore ^ -rr = 2 -j- 17 

a dt d8 dt 

where — Vr- denotes the differential coefficient of R with 
dt 
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respect to t, only so far as It involves t through involving the 
co-ordinates and elements of the disturbed planet. Now since 
the velocity of the planet at any time can he expressed in 
terms of the co-ordinates and elements of the instantaneous 
ellipse constructed for that time in the same form as if it 
moved in this ellipse, its component in any direction can also 
be so expressed. Hence, considering ii as a function of 

df du dz 
r^, 0^, z (see Art. 11), the values of -ttS ~^ ^'^^ ^^'^ there- 
fore of \, , may be expressed in the same forms as if the 

elements were invariable. Since, then, t always occurs in R 
coupled with e in the expression ni + e, we have 

d (B) dR dR 



therefore 



dt d {nt 4r e) de * 

da 2na^ dR 



for 



dt fJL de * 

27. To obtain aformtdafor calculating the excmtricity. 
From equations (3) and (4), if w^, be eliminated, we obtain 

\dt) " h' [r h) • 
Diflferentiating as if r were constant, and writing --r- 

de' 

dh 
dt 



smce 



dr 
dt' 


fiede 
~ h' dt 


{fHi-f.)C-?)) 


s 


fi'ede 
~ h* dr 


h* 


■ e')dB hdR 

dd t* de' 


dh 


dR 




9 


dt' 


~ dd' 
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, f ii^e ds _dR dr dR de^ fi^l-e^dR 
tneretore ^. ^^ ^ __ + _ _ ^,_ _ 

fdRdr dR d0\ //>' (1 - e*) dJi 

.1. /• cte nh*dR 1 — ^/dR . c?B\ ... --v 

therefore ^«__ - -_(^_ +_j, (Art. 16), 

/(ie de fie \de dw)^ 

since A' = /*a (1 — c'), and nV = /Lt. 

28. This formula may also be deduced from that of the 
mean distance, by means of the equation 

A' = /xa (1 -- e*). 
Wehave 2Ag-,.(l ^6«)g-2/.a6 g, 

or 2naV(l-e^)^= 2na« (1 - e») ^ - ^/ioe^^; 



29. To obtain a formula for calculating the longitude of 
perihelion. 

From equations (3) and (4), if e be eliminated, we obtain 

dr . f^ . h fJL 

Dififerentiating as if r and 6^ were constant, and writing 

dR r dfr , . . 

-^ for -Tsr , we obtam 
dr dt 



^cosec'(^,-.t^,)-^Hcot(tf,^i^,)^-(^- + pj^ 

1 /xXdfi 
r'^hVde ' 
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therefore ^ cosec (0^ - tt J -^ 



dR , /I , ti\ . .^ .dR 

dd 



= -cos(^,^^,)g + (J + gsin(^,-t.,) 



1 dR / A , , h*\ 

but from equation (4) ^ cosec (^^ — ^J"©) ~ ^ > 

, « ^^o_ ^ dR ^natJ(l — e^dR 

dt fiea de ^ fie de * 

Now if tir denote the longitude of perihelion measured on 
the plane of reference as far as the node, and thence on the 
plane of the orbit, fl the longitude of the node on the plane 
of reference, Hg its longitude on that of the orbit, we have 

«r — -oTj^ = fi — il^; 

therefore d^ ^d^, ^dil dn, 

dt dt dt dt * 

Now -17 is the angular velocity of the line of nodes on 

the plane of reference, —^ its angular velocity on the plane 

of the orbit ; and since the plane of reference is fixed, the 
former is the total angular velocity : hence 

da. dH . 
W^lt''''> 

therefore -3- = -yr -f (1 — cos i) -j- : 

dt dt ^ ^ dt 

or, substituting for -1— from Art. 31 or 36, 

, <far^ r»oV(l-e') dE ""*^2 dig 
dt fie de fi V(l — «*) di ' 
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To obtain formulas for calculating the longitude of the 
node, and the inclination. 

30. We now return to our third equation of motion, 

. dO^ » 

or, as it may be written, 

% = Zr (1). 

We have seen (Art. 18), that <^a = 0; hence the motion of 
the plane of xy, which coincides with the plane of the orbit, 
is compounded of the angular velocities ify^ about the axis 
of X, and ^3 about the axis of z. Now the former is equiva- 
lent to an angular velocity (f)^ cos {0 — O) about the line of 
nodes, and an angular velocity <^j sin (0 — fl) about an axis 
perpendicular to it in the plane of the orbit : but the angular 

di 
velocities of the plane of the orbit about these axes are -r and 

sin i -^ respectively ; therefore 

(f>^ cos (^ — ^) = ;77 , ^1 sin (d — fl) = sin i -j- • 
Hence, by equation (1), 

Ag = Zrcos(^-fl) (2), 

Asini^ = Zrsin(d-fl) (3). 

31. In order to determine Z we must suppose the curve 
of reference perpendicular to the plane of the orbit. If we 
denote by 5 an arc of this curve measured from some fixed 

d R 
point up to the planet, we have by Art. 8, Z= -p • Now 

in the same way that the position of the planet is known 
when we know r, 6, i and il, the position of any point on 
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the curve of reference may be determined by its polar co- 
ordinates r, on a, plane passing through it and the Sun, the 
inclination i of this plane to the plane of reference, and the 
longitude fl of its node. Since, however, an infinite number 
of planes can be drawn through two given points, we must 
introduce some further condition to fix the position of that 
on which r and are measured*. 

Dififerent forms of expression will be obtained for Z 
according as dififerent conditions are assigned. First sup- 
pose the plane to pass through SN, the planet's line of 
nodes : let P be the position of the planet, draw FG per- 
pendicular to SN, and take for the curve of reference a 
circle AP with centre C. Through SN draw a plane in- 
clined at a small angle &i to the plane of the orbit, cutting 
the circle AP in p, and let Pp==B8; then since Pp is per- 
pendicular to the plane SNP, 

Pp^CP.Bi 

Bs^r sin (0 — fl) Si ; 



or 



therefore 



di 



da r sin {0 — ft) 




See Appendix, Art. 7, 
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Now, (see Art. 13), B may be expressed in the form 

i?=/Cr,^,n,i), 

in which, if the series for r and be substituted, B will be 
expressed as a function of t and the elements. Hence 

ds " dr ds dO ds dil ds di ds * 
, dr ^ d0 ^ dSl ^ di ^ 



therefore 



ds ' ds ' ds ' ds rsin(tf-fl)' 

dR 1 dR 

ds ~ r sin (5 — X2) di ' 



On substituting this value for Z in equation (3), we 

obtain 

, . .dn dR 
ti sm % —T7 = -^n 1 
at di 

or since A* = fta (1 — e'), and nV = /i, 

rf^ /ivll^O^iii* ^* 

32. By substituting the value of Z found in the pre- 
ceding Article in equation (2) of Art. 30, we might, of 

di 
course, obtain an expression for -^ , but this would involve 0, 

and consequently be in a form inconvenient for calculation. 
We proceed^ then, to obtain an expression for Z by means of 
which may be eliminated. 

33. Suppose the plane on which r and are measured, 
instead of passing through SJS', to pass through a line SC 
in the plane of the orbit perpendicular to SN : let P be the 
position of the planet, draw PC perpendicular to SC, and 
take for the curve of reference a circle AP with centre C. 
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Through SC draw a plane inclined at a small angle to the 
plane of the orbit, cutting the sphere in the great circle nmp, 
and the circle AP in p. Draw Nm perpendicular to np. 




Then Nm which measures the inclination of the two 
planes == ^ SO sin t, and 

CP=rcos(5~fl); 
hence if Pp = 8$, 

& r= — r cos (^ — fl) Sft sin i; 
Sil -1 



therefore 



Again 



therefore 



also 



Hence 



dR 
ds 



ds r cos (0 — fl) sin i ' 

S0=^Ln + np-'{LN+J!i'P) 

= nm — nN" 

= — Sfl cos t + Sfl ; 

d0 ,^ . dil 

^ = (l^cost)^; 

- = ^ = 0. 

dBde dRda 
dd ds da ds 

[dR 



[dR , ,- ., diJ) 



r cos (tf — X2) sin i 
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On substituting this value for Z in equation (2), we have 

1 fdJB . « . ^i fdR , dR\] ,.^ _, 
— Arol^ + 2Bm«-(^+^)} (Art. 15), 

na f 1 dR i(dR dR\] 

The only remaining element is the epoch, but before pro- 
ceeding to obtain a formula for its calculation, we shall give 
another method of obtaining the results of Arts. 31 and 33. 

34. To obtain a formula Jbr calculating the inclination, 
{Second method,) 

If the motion of the planet be referred to the polar co- 
ordinates of its projection on the fixed plane of reference and 
its distance from this plane, we have the equation 

1 d f ^de^ IdR .^ . . Q. 
7JtV'^'di)^r,de,^^''^'^'^^^ 



or 



rf/ ^d0A_dR 
dt\^ dt) d0,' 



Now if SA, BA^ denote the vectorial areas swept out in 
the time Bt on the plane of the orbit and the plane of reference 
respectively, we have 



BA^ = SA cost: 



but BA,^lr,'B0,,SA=\r'Bd,; 

therefore r," -37 = ^* -jt cos i = h cos i. 

^ at at 
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or 



Hence our equation of motion becomes 

d ,. .. dB, 

, . »di ^ ,dh dR 



^y t. 1 . »d% dR ,dh 

therefore — A sm i j- = -37^ — cos i -f- 

dt da^ dt 



dR 



»dR 



^mr'^'dd 



therefore 37 = 
db 



dR , ,, ..fdR , dR\ ,.^ ^^, 

=._ + (!, eos^)(-^ + ^j, (Art. 15); 

1 (dR . ,, . fdR . dR\] 

~7r~^ — = 1:70 + (1 -cos*) rij- +T- r 
A sin t [dil ^ ^ xde d'srj) 



na 



( 1 dR 



fi V(l — ^) (sin i dD. 



i (dR . dR\\ 
dm)] 



+ *^^2(,^+--'^- 



35. To obtain a formula Jhr calculating the longitude of 
the node. (Second method.) 



We have by Art. 13, 

ii=f(r.e,n,t) 

or, since evidently — Q,='0^—il^, 



(1). 



(2), 



tiff 

* We have made this transformation^ because, although the value of -j^ 

at 

is the same in form as if the elements were invariable, this is not the case 

de 



with 



dt' 
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Now we have seen (Art. 26) that \ (denoting by it 

that R is to be diflferentiated with respect to t only so far 
as it involves t through involving the position of the dis- 
turbed planet) may be expressed in the same form as if the 
elements were invariable. 

We have then, considering the elements variable, 

d{R) ^dRdr c?Jg d (g^ + Q - Qp) dRd£i dRdi 
dt "drdt'^dd^ di "^ da dt "^ di dt* 

and, considering them invariable, 

d(R) ^dRdr dRdO^ 
dt dr dt dO^ dt ' 

Equating the two values of ^ , we obtain 

dR (do. _ dflA .^^,^^'^A 
doXdt dt J "*"dn dt "^ di dt 

-- dR . /a\ ^^ • /t\ dR . dR 

Now ^^m(2)=-^m(l)=;^ + ^. 

and (see Art. 29), -^ = ^ cos i ; 

therefore |^ + (1-C08,) (^ +^j}-^ +^ ^ = 0. 

di 
Substituting for -^ its value, 

dj£l _ na dR 

dt " fi V(l — «*) sin i di ' 

36. To obtain a formula for calculating the longitude of 
ike epoch. 

If R be expressed in terms of t and the elements (see 

Art. 13), since w* + e always occurs as one symbol, we may 

write 

R =f{nt + €, a, e, «r, il, %). 



■ 

* 

I 
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Diflferentiating, the elements being considered variable^ 
we have 

d{R) dRd(nt + €) dRda dRde dRdfsr 
dt de dt da dt de dt dtsr dt 

dRd£l dRdi 
■^dXl dt "*" di dt' 

and differentiating as if the elements were invariable, which 
is permissible for the reason explained in Art. 26, 

d{R) ^ dR 
dt de * 

Equating the two values of , , 

dR _ dR / dn d€\ dR da dR de dRdm 

de d€\ dt dt) da dt deHi dtsr dt 

dR da dR di 
dH dt di dt' 

Substituting for -^ , -^ , &c., their values 

0= — /'* — + — ^ ^^' ^^ ^^ rta (1 - g*) dR dR 
de \ dt dt) /i de da fie de de 

nas/il-^) [dR , dR\dR , na^/ii-e") dR dR 
fie \d€ dwj de fie de dm 

"°^^2 dRdR na dRdR 

,. ./*-y'(l — e') di dvr fj,»/{l — e'ja.ni di dil 

9 

__ na dRdR ^^^2 /dR dR\dR . 

fi*J{l—e^)WLidQidi fif^(l—e^)\de dur) di ' 

c. p. T. 3 
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^dRf dn d€\ 2»ia' dR dB 
de \ dt dt) II de da 

na ^/(l — 6*) ,- , .- «x , dR dR 

"«**^l dRdR 



fA V(l —ej de di ' 

Dividing every term by -r- , and transposing, we obtain 
de dn inc? dR naV(l— c*),- ,„ ^^^dR 

37. Of the formulae which have been obtained for cal- 
culating the elements of the orbit, that of the preceding Article 
is the only one which contains a term proportional to the 
time*. It may, however, be replaced by one in which no 
such term exists. For, let f denote the mean longitude, 

then 

f = 7i< + € ; 

therefore | = „ + ,| + |, 

^ /^ r jj\ .dn ^ de 
or ^,(S-!rdt)^t^^^j,^ 

Now let f=/nd^ + e', 

then by the formula of the last Article 

de 2na^ dR . na ^/(l - O r^ ,., _«., dR 

«atan£^ 

'^(i>/{l-^ di ' 

* The reader, if acqnainted vith the Lnnur Theoiyi vriU have already seen 
the inoanTenienoe of snoh terms. 
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Since in the elliptic formulse € never occurs except when 
coupled with nt, in the expression nt + €, it will be altogether 
eliminated if for w^ + € we write fndt + /. 

Considered as replacing the element €, e' is called the 
epoch*. Since however we shall never have occasion to em- 
ploy the formula of the preceding Article, the accent will in 
future be omitted. 

fndt is termed the mean motion in the disturbed orbit, 
and is denoted by ^, 

38. To obtain a formula for calculating ^, we have 

^^dn 
df dt' 

and diflferentiating the equation nV ^ fi, 

dt dt 

- « dn^ Sn da _ Sn*a dR 

dt" 2a dt" fjk de* 

cfC^ 3n'a dR 
d^ ji de ' 

39. We will here recapitulate the formulsB which have 
been obtained for calculating the elements of the orbit : 

^.» da _ 2na* dR 

de_ na{l-^) dR na^/jl"^) (dR dR\ 
^^^^ dt tie de ps Ue rfW' 

..... du _naWO^-f) dR J^^^^^dR 
^*"^ d« ~ /JUS de ti'/(l- e") di ' 

* It may be notioed that if n were oonstant, e' woold be identical with e. 

3-2 
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(V) 



natan| ^ 
dD, na dR 



dt fi V(l — O sin i di ' 

^ .. di_ na ( 1 dR i fdR dR\\ 

^""'^ de" /i V(l -«") tsini dXl ■^'^^'^ 2 Vd^ "^ dWj ' 

We have also (vii) the equation 

d'C^ 3n*a diZ 

but this forms no new relation, since it has been deduced 
from (i). 

40. When the elements have been calculated by means 
of the above formulae, the position of the planet will be given 
by the equations 

r = a]l +^e' — ecos(5'+€ — 'Br) — ^e'cos2($'+6 — -cj)— ...^ , 

5 

^ = (;+6 + 2esin(C+e-«r) + 7e"sin2(^+e-'Br) + .... 

NoUi, Some of the Author's methods igre similar to those given by 
Mr M. O'Brien in an article **0n certain formula in Physical Astronomy,** 
published in the Camb. Math, Journal^ 1843, Vol. iii. pp. 249 — 259. Com- 
pare sections 18, 29, 30, 8, 20, and 31 of this work with sections 6, 7, 8, 10, 
11 of the article releired to. [a. f.] 



CHAPTER III 

DEVELOPMENT OF THE DISTUKBING FUNCTION. 

41. In the first Chapter we have obtained equations by 
means of which B may be expressed in terms of the time 
and the elements of the orbit ; we now proceed to shew how 
the actual development may be effected in a series ascending 
by powers and products of the excentricities and the tangents 
of the inclinations. In the Planetary Theory these are ex- 
tremely smaU, and the series wiU converge rapidly. Ac 
cordingly in the present treatise small quantities of orders 
liigher than the second will be neglected*. 

42. If we recur to Art. 11, it will be seen that, consider- 
ing only one disturbing planet, 

iJ = mT ^ r 

L{r,« + r;« - 2r,r; cos (6^ - O;) + (^ - /)»}* 



r^r/ cos (g,-g/) + ^^n 



The first step towards the required development will be 
the expansion of r^, r/, d^y 6^, z and / in terms of the time 

* We may remark that to this order of approximation the inclinations, 
their sines, and tangents will be equal 
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and the elements of the orbit. For this purpose we may 
employ the equations which have already been obtained in 
Art. 13, viz. : 

tan {0^ — il) = cos t tan (0 — 11), 
sin \ = sin i sin {6 — il), 
r^ = r cos \, z^r sin \, 

r=sa-|l + -^e* -ecoa(nt + €- w) --^e^ cos2(nt + € — vf)'-...i , 

5 

^ = n^ + € + 26sin(n^ + € — 'Bj) + T«*sin2(n< + € — w) + ..., 

with similar equations involving the co-ordinates and elements 
of the disturbing planet. 

(i) To expand r^. We have 

r J = r cos \ = r (1 — sin' X)* 

ar f 1— 5 8in*\ + ,..j 

= r Jl- itan^•sin»(^-n)+...l 
to the same order of approximation, 

f 1 1 . 

= r -^1 — T tan* i +7 tan' t cos 2 (5 — il) + ...^ ; 

or substituting the expansions for r and 0, 
f 1 1 . 1 

+ ~ tan' i cos 2 (w^ + € — fl) + . . . !■ 
= a (1 + m), suppose. 

Similarly, r/ = a (1 + %')• 



==^••11— ^sin'tsin' 
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(ii) To expand 0^. We have 

tan(d,-^) = tan{(^,-n)-(^-n)} 

tan(g^--ft)-tan(g~fl) 
1 + taii (^j - XI) tan (^ -O) 

_ (cos i — 1) tan (0 — fl) 
"" l + cositan"(^ — ft) 

• 
-2siii''|taii(^-li) 



1 + tan* (^ - II) - 2 sin* | tan'(^ - Q) 



-sin*|sin2(^-n) 



l-2sin*|8m*(^-fi) 



»-8in*|8in2(^-li)-...; 
therefore ^^ — ^ = — sin* ^ sin 2 (^ — ft) — ... 

= - tan*|sin2 (^- ft) - ... , 

to the same order of approximation ; or^ substituting the 
expansion for 0, 

^j =s «< + € + 2e sin (n^ + € — tj) + T 6* sin 2 (n< + c — «r) 



— tan* ^ sin 2 (n< + € — ft) + . 
^nt + € + v, suppose. 
Similarly, ^/ = n't + e' + v\ 
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(iii) To expand z. We have 

-y = r sin X = r sin i sin (5 — ft) 

= r tant sin(5 — ft)— ... 

to the second order ; or, substituting the expansions for r 
and 0y 

z = a {tan i sin (n^ + e — ft) + ...], 

in which terms of the second order will not be required. 

A similar expression may be found for /. 

43. Having obtained the expansions of r^, r/, 0^, 6^\ ^ 
z, / we must now substitute them in. the expression for JB. 
This may be effected as follows. 

Let JB, be the value of R when u, u, v, v' are severally 
zero : then, writing ^ for nt +€ — (n't + e') we have 

B^ = m' [{a» + a'* - 2aa cos ^ + (^ - /)»}-*. 

- {aa cos ^ + zz') (a"^ + z"^"^] 

= m [(a* + a'* — 2aa' cos (j>)~^ — ^ cos if>] 



- w' ll (a* + a'* - 2ao' cos ^)-» (^ - /)* 



^ 1 , 3 a ^ .1 



N w JJ = -B^ + -^ au + ^/ a V + ^ (^ - « ) 



+ 









+ 
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44. It will be shewn in a subsequent Article that. 

(a* + a" — 2aa cos ^)"' 
can be expanded isL a series of the form 

o .^, + jIj cos ^ + A^coa2^ + ...+ Aiicoa k^ + ... 

Assume then 

(a* + a'*-2aa' cos <f>)~i = „ C, + C, cos ^ + C, cos 2^+ ... 

(a* + a" - 2aa' cos ^) " * = 5 -Dj+ Dj cos ^ + 2?, cos 2^ + . . . 
Thus 
B = m' II Co + f(7, - ^1 cos ^ + (7, cos 2^ + ...I 

^ , , , fl do, , fdC, 2a\ ' , d(7, „ . . 1 
-m{v- v') |(c7j - p) sin ^ + 2(7, sin 2^ + . ..I 

__ m-(i>-t,T |^Cr,-^,)cos^ + 4C,co3 2^+...] 
-m'o«(t;-t;7K^-l)8in^ + 2^sin2^+..j 



42 
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-«.'aV(«-t,'){(g + ^)siu^ + 2g»sm2^4-...} 



m!{z-£y fl 



J5 2)o+ DjCOS^ +2),C08 2(^ + ...j' 






45. By Art. 42, 



M = ~e* — -7 tan*t — 



1^ 



e COS (w^ + € — 'ct) — ^ e* COS 2 (n^ + 6 — «r) 



+ ^ tan' i cos 2 (n^ + e — fl) + . . , 



5 

t; = 2e sin (n^ + € — «r) + j e*sin 2 (nf + e — -or) 

— tan'^ sin 2 («^ + € — fl) + ..., 
J? = a { tan i sin (w^ + € — fl) + . . . } , 
with similar expressions for t^', v', /. 

Hence w' = e*cos* (n^ + € — -or) + ... 

(v - 1;')« = 4e« sin* (w* + e - tcr) + 46" sin" (n'« + e' - «r') 

— See' sin (w^ + e- -cr) sin (n ^ + e' — «r') + . . . 

= 2 (e* + O - 2e* cos 2 (n^ + € -- tsr) - 2e'" cos 2 {nt + e' - -or') 

— 4ee'cos(<^-'Br+«r') + 4ee'cos {(n + n') ^ + e+ e'— «r — 'or'}+..., 

W = ee'cos (w^ + € — -cr) cos (w'^ + c' — -cr') + ... 

e^ ee 

= -^ cos (<^ - «r + «r') + -^ cos {(w + w') * + € + € - w — «■'} + . . . 
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u (v — 1;')= "" ^ sin 2 (n^ + e — 'cr) 

+ 2^e' cos {nt + € — -cr) sin (n't + c' — 'cr') + ... 
= — e' sin 2 (n^ + 6 — «r) — ee' sin (^ — -or + «r') 

+ ee' sin {(w + n') ^ + c + e' — -cr — -or'} + ... , 

(«--8^) = — 2 — ■*■ — 2 2 — cos2(w^ + 6-fl) 

a'*tan'i' 
^ — cos2 (nt + e— OTf—aa tanttan i'cos(<^ — fi + Of) 

+ aa! tan /tan i' cos {(n + n') < 4- e + e' — fl — fl'} +... 
&c. 

46. If these values be substituted in Art. 44, it will be 
seen that cosines will be multiplied only by cosines, and sines 
by sines. Hence the series will consist of two parts, one inde- 
pendent of t explicitly, and the other consisting of periodical 
terms of the form 

P co^ {{pn±qn')t+Q], 

where p and q are any positive integers or zero, P is a function 
of the mean distances, excentricities, and inclinations, and Q 
a function of the longitudes of perihelia, nodes, and epochs. 
The former part is denoted by the symbol F: we proceed to 
determine its value as far as the second order of small quan- 
tities. 

47. To determine that part of R which is independent of 
the time esopliciily. 

If those terms only be written down which either are, or 
after reduction will become, independent of t, we have 

, f C, , g dC^ (^ tan' A a' dG^ /e^ tan't\ 
^ (2 "^2 rfa V2 4 r^ da'V2 4 / 

a'd^o^ a:'d?C,e'^ D ja^\AuH a'tanV X 
"*■ 4 (fo» 2 "^ 4 da"' 2 4 V 2 ^ 2 J 



m 
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+ 2 f 0^ Taj COS (f> 4iee' cos (^ — -or + «r') 



+ a f -^ Tsj sin ^ ee' sin (^ — «r + «•') 

^ l^ "*" ~^J ^^^ 9 ee sin (^ — «r + «r ) 
+ ^ cos ^ aa' tan i tan t ' cos (^ — fl + ft') + ...[ . 

Now cos ^ cos (^ — «r 4- «r') and 8in0sin(^--'CT+«r') con- 
tain the term J cos (isr — -or'), cos (f> cos (^ — fl + fl') contains 
the term J cos (fl — fl') ; hence 

+ 2 cw'Di tan t tan i' cos (fl — fl') + ...V . 

We shall hereafter be able to simplify this expression. 

48. We have seen that the remaining terms of It are 
of the form P cos [(pn ±qn') t + Q]: if then values of p and 
q could be found such that ^ + qn' = 0, this term, being in- 
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dependent of t explicitly, would form an additional term in F. 
No instance of this, however, occurs among the planets. 

49. In consequence of the extreme smallness of the ex- 
centricities and inclinations of the orbits of the principal 
planets, terms in B of orders higher than the second may in 
general be neglected : but it sometimes happens, as in the 
Lunar Theory, that higher terms become sensible through 
the process of integration. This we shall consider in a sub- 
sequent Chapter, but the following proposition has an impor- 
tant bearing on the subject. 

50. The principal part of the coefficient of a term in R 
of the form P cos {(pn — qn') t + Q} is of the order p ~ q. 

Def. By the principal part of the coefficient is meant 
that part of P which is of lowest dimensions in e, e\ tan t, 
tan i\ 

If we return to the expression for B, in Art. 44, it will 
be seen that in order to obtain the general term it will be 
necessary to multiply the product of the general terms of the 
expansions for i^*, w'^, v^, v'*, ^, z^ by cos Jc^ or sin fc^. 

Now (i) in the expansions of w, w', t;, v', z, z the follow- 
ing law is observed to hold : — The number which multiplies 
nt'\'€ or nt + e in the argument of any term represents the 
order of the principal part of the coefficient of that term. 

(2) The same holds good in any power of w, u\ v, v\ 
2f, or z\ Tor consider a term Pcos (pnt-i-q) in u\ It can 
only have arisen in the following ways; partly from the 
multiplication of two terms in u of which the arguments are 
Int + \ and mnt + fi, where l + m=p; and partly from such as 
have the arguments Vnt + \' and m'nt-\-fi, where Z'~m'=^. 
In the former case the order of the coefficient will be l + m, 
which equals p, in the latter it will be V + m\ and this is 
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greater than p. Hence the principal part of the coefficient 
of a term P cos {pnt + q) in u*, will be of the order pi 

Since then the law holds in u^, it may be shewn in like 
manner to hold in the product of w* and u, i.e. in u\ Thus 
it may be proved for any power of u. In like manner it 
may be shewn to hold for any powers of u, v, v\ z or z\ 

(3) The same law is true for the product of any powers 
of M, v, ^; and likewise for the product of any powers of 
«', v\ z\ This may be proved by a method similar to that 
of (2). 

(4) In the product of any powers of w, w', t;, xi\ z and 
z\ the order of the principal part of the coefficient is the 
arithmetical sum of the multipliers of nt and v!t. 

For let us consider a term Jf cos {(Zw ± ZV) t + N\. Now 
this must evidently have arisen from the multiplication of 
L cos {pnt + X) with L' cos (ZV^ + X')> or of L sin i^ + \) 
with L sin {Knt + X'), where by (3) L is of the order I and 
L' of the order l\ Hence M will be of the order I + 1\ 

Now any term in the development of R of the form 

Pcos {(jpw — j«') < + ^} must have arisen partly from the 

cos 
multiplication of P^ . fc^, or as it may be written 

with P, '2. {[(p _ ^) „_ (j _&)«']< + Q.}. 

and partly from its multiplication with 

where h is any positive integer or zero, P^ is a function of a 
and a' only, and P^, P3 are functions of the excentricities and 
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inclinations, such that the orders of their principal parts iEure 
given by law (4). Hence the order of the principal part of 
P will be equal to the lesser of those of P^ and P^. 

Now the order of the principal part of P^ will be the 
least value of which the arithmetical sum oi p-^k and j - i 
is susceptible, for different values of k. 

(i) Suppose k intermediate to p and q\ then this sum 

= J? - & + g - A: =^ - g ; 

(ii) Suppose k not greater than the smaller of p and q ; 
then this sum ^p + g — 2A?, the least value of which (by 
putting k equal to the smaller of jp and q) = j? - q ; 

(iii) Suppose k not less than the greater of p and q\ 
then this sum= 2Z;— 1> — y, the least value of which (by 
putting k equal to the greater oip and q)=^p'^ q. 

Thus p - g is the order of the principal part of P,. Thai 
of Pj will be the least value of which y + fc + g + A is sus- 
ceptible, i.e. p + g. 

Hence it appears that the order of the principal part of 
P \ap^q. 

51. The principal part of the coefficient of a term in R, 
of the form P cos {(pn + qn')t + Q} is of the order p + q. 

This term arises from the multiplication of such terms as 

with P.'Xj.{v-'k)n^{l^'k)n'\t+Qi,, 

and P, "2, {[(i* + *)« + (?- h)n^t + Q.}, 
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and as in the last Article, the order of the principal part 
of P will be equal to the lesser of those of P, and P,. 

Now the order of the principal part of P, will be the 
least value which the arithmetical sum oi p^^h and q-\-h 
can assume, for different values of h 

(i) Suppose h less than p\ then this sum 

= J? — i + J + i = J9 + J. 

(ii) Suppose h not less than p ; then this sum 

the least value of which (by putting k equal to p) =p + gr. 

Similarly it may be shewn that p + q will be the order 
of the principal part of Pg. 

Hence it follows that p + q will be the order of the 
principal part of P.. 

In Art. 44 we have assumed that (a* + a" - 2aa cos ^)"* 
can be expanded in a series of cosines of (f) and its multiples, 
we shall now give a proof of this and shew how the coeflS- 
cients may be calculated. 

52. To shew that (a' + a'* — 2aa' cos ^)"* can be expanded 
in a series of cosines of multiples of<f>. 

Suppose a greater than a', and for - write a; then 

a 

(o" + a" - 2aa'<os <!>)'• = a^{l + a.*- 2a cos ^pp 
= a-* {1 + a» - a (c*^~^ + e-*^}-' 
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d^il+82e + ::: — ^ a V 



2 



W.9 + l)(g + 2 ) .j^^-i . 

+ Tn — " ^'^ + • • • 



X < 



1 + s:.e + - — ^ — 0.6 



+ l(i±i)J-l±-?la»e"'*^'^ + ... 



= a-''ri+sV 



+ 



a 



2 



+ 



[s(s + l)(s + 2) 



15 



JV+...1 



+ 2a^^sa+s.?^j"-)a' 



i 



+ 



s(5+l) g(^+l)(g+2) 
2 • 



a^+... 




e +e 



) 



+ 



where the coefficients of e and e will always be 

equal. Hence we may write 

(a'+a''-2aa'cos^)"* 



where -4^, A^, &c., are functions of a and a. The series 
which they represent will be always convergent provided a is 
less than unity, or a greater than a\ If a be less than a, we 
have only to interchange a and a' in the above, so that a will 
then denote the ratio of a to a\ 

c. P. T» 4 
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53. To calculate G^ and C^. 

In the preceding Article, let 5 = ^ ; then 

^ „_, (1 11. 3, ,1.31. 3. 5,, ) 

C,= 2a |-a + ^.2-^a+2-^.2--^ga+...J. 

Unless a be small, these series will converge too slowly 
to be practically useful. More convergent series might be 
obtained, but Pont^coulant remarks (Systeme du Monde, 
Tome III. p. 81), that in practice it is more convenient to 
employ elliptic integrals for the purpose, in the manner we 
proceed to explain. We have 

(a* + a ' - 2aa cos ^y ^ 

^^ i = |acos6+i(7,(l + cos2<^) 

(a' + a'»-2aa'cos<^)* ^ ' ^ 2 '^ ^^ 

-f- X C7j (cos 3^ + cos ^) + . . . 

Integrating both sides of these equations with respect to <f> 
between the limits and 2*rr, we obtain 



TT 



' Jo {a'+a''-2aaco3<f>)^ «^o (1 +a»~2acos^)*^ ' 



^ __ r2» cos <^ d(l> _ 1 P" cos (f> d<j) 

'"Jo (a« + a"-2aacos^)*"«Jo (1 + a' -2a cos <^)*' 

These integrals may be reduced to the standard forms of 
elliptic functions by assuming 

sin(^ — ^)=asin^ (1), 

whence tan^ = — . ^ (2). 

CQS 9 — a 
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From (1) cos {6 - ^) (l - ^|) = a cos 0) 

therefore ^ - cos (g - <^) - a cos g 
theretore ^^ - cos(g-(^) * 

Now 
cos (g — ^) — a cos g = cos 6 (cos ^ — a) + sin fl sin ^ 

/cos' g . . /,N . . , /ON 

sin <^ 
sin 6^ 

= V{(cos^-a)« + sin'<^l, by (2), 

= V(l+a'-2acos^) (3). 

Also cos (5 - ^) = V(l - a' sin'* ^) (4). 

TT d6 //I 4- a' - 2a cos 6\ 

^^^"" ^ = V l 1-a-sin-g J' 

Again, from equations (3) and (4) 

V(l + a»- 2acos <f>) = V(l - a' sin' ff) -a cos5; 

therefore 1 + a' — 2a cos ^ = 1 — a' sin' 6 

+ a' cos' ^ - 2a cos 6 *^{1 - a' sin' 6), 
2a cos <^ = 2a' sin' ^ + 2a cos V(l - a' sin' 6), 

or cos ^ = a sin' + cos ^ \/(l — ot* sin' ^). 

Now as <^ increases from up to 27r, ^ also increases from 
to 27r ; hence 

" aTTJo V(l + a* - 2a cos ^) 



OTTJo V(l-a*8ia'6')' 

4—2 
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^ _ 1 P'^ COS j> d(f) 

^ air Jo V(l + ci* — 2a cos <^) 

1 P' asin'^dd . 1 /*' /, ,/, 
= — / Tr-t — 8 • « /i\ H / cos tfau 

=— ir-/ /i ?» 8^x -rva-ct»siD«g)d4> 

a27r yo V(l-» sm'^) Jo J 

Hence with the usual notation for elliptic integrals (see 
Todhunter's Integral Calculus, Art. 222), 

The numerical values of F(a,-^] and j&fa, ^j may be 
found from Legendre's tables of elliptic functions. 

54. CUven C,, and G^_^ to obtain C,,^^. 
We have 

(a^+a"^^ 2aa'cos ^)"*= « ^o"*" C^iC<^8^ + ••• + C^»cosA?^ + ...; 

dififerentiating with respect io <f>, 

aa sin ^ (a* + a" — 2aa' cos ^)"^ = C^ sin ^ + 2 (7, sin 2^ + . . . 

+ A;Cj^sinfc^ + ...; 

therefore ad sin ^ ( ^ C^ + (7^ cos ^ + . . . J 

= {g? + a*- 2aa' cos ^) (C^ sin <^ + 2(7,sin 2<^ + ...); 
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equating coefficients of sin k<l>, 

I aa'iG^-C^,)=k{a'+ a") C,-aa'{{k - 1) G^, + (fc + 1) C^,}; 

, ^ 2fc a^ + a'^ ^ 2k-l ^ 

whence w.i = of -i 7 — G^ — t,i — -i C^_,. 

"-^ 2*4- 1 aa * 2A + 1 *"* 

55. (riven C^ and 0^+1 to obtain D,,. 

As in the last Article, we have 

aa' sin (f) (a* + a' — 2aa' cos ^)'* = (7^ sin ^ + 2 (7, sin 2^ + . . . 

+ fcCjSinA;^ + ... ; 

therefore aa' sin ^ ( ^ D^ + D^ cos ^ + D, cos 2^ + . . . j 

= (7jSin^ + 2C4sin2^ + ... ; 
equating coefficients of sin kif), 

2A;(7,aa'(A.x-^*+i) (l)' 

writing fc + 1 for k, 

2(Je+\)C^,=^aa'{D,-D^) (2). 

Again, 

(a' + a'*- 2aa' cos <^) * = ^ D^ + A ^^^s ^ + A ^^s 2<^ + . . . , 

and 

(«• + a"- 2aa' cos <^)"* = ^ ^o+ ^i^^s <^ + (7^ cos 2^ + ... ; 

therefore ^ 0<, + (7^ cos ^ + . . . 

= («' + »''- 2aa'cos if>) ^gD^ + J^^cos <^ + ...j : 

equating coefficients of Cos fc^, 

0,= (a* + o'«) A-aa'(Z)^, + A„) (3). 

writing A? + 1 for i, 

Cr^. = (a' + a'')2)^,-aa'(A + 2)»J .......... (4). 



54 PLANETARY THEOEY. 

Eliminating 2)j_j between (1) and (3), 

{2k+l) C, =; (a' + a'^D,- 2aa'i>^, ......... (5). 

Eliminating Dj^^ between (2) and (4), 

(2*4-1) a^, = -(a» + a'')A,,+ 2aa'i>, (6). 

Finally, eliminating Dj^^ between (5) and (6), 

■ 

(2k + 1) {(a* + a'») C, - 2aa' C^,} = {(a* + ay - 4a*a''] D„ 



or A = 



(^a-Jay K«* + a'O C', - Zoo'C^J. 



56. To calculate the successive differential coefficients of 
C,, and Jy^ with respect to a and af. 

We have 

differentiating with respect to a, 

— {a-' a cos (j)) {a? + a'* — 2aa' cos ^) = ^ -r-** + -i-^ cos^f . . . 

+ -T-* 008*9+ '"J 

substituting for (a* -f a'* — 2aa' cos ^)~* its expression in series 
— (a — a'cos^) f^2>^+DjCos^+ ... +i>jCos*^+ .,. j 

1(^(7, dC^ JL . . dC^* 7JL . 
= 7:-,-*^ + -7-^ cos 6+,., +-5-* cos &<!>+...: 
2 da da ^ da, ^ 

equating coefficients of cos h<f>, 
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By giving to A; in succession the values 1, 2, 3, &e., those 

dC dC 
of --j-^ , -j^ , &c. may be found, the right-hand member 

being calculated by the formula of the preceding Article. By 
equating the parts independent of ^, we obtain 

The value of -j- may be found by dififerentiating the 

dC 
expression for D^ in Art. 55, and substituting for -^ and 

dC 

-^ their values as given by the present Article. 

The successive differential coefficients of (7^ and Dj^ with 

dC 
respect to a may be obtained from the expressions for -7-* 

and ,-^ by simple differentiation and substitution. 

57. We might determine in the same way the successive 
differential coefficients of C|^ and Dj^ with respect to a' ; but 
when those with respect to a have been found, the former 
may be derived from them, as we proceed to shew. On ex- 
amining the expansion of (a* + a!^ — 2aa' cos 0)"* in Art. 52, 
it will be seen that Aj^ is a homogeneous function of a and a' 
of — 2^ dimensions. Hence C^ and Z)^ are homogeneous 
functions of a and a', the former of —1, the latter of —3 

dimensions. It follows that -r^ , -j^ will be homogeneous 

functions of — 2 and — 4 dimensions respectively ; and so on. 
Now by a known property of such functions 

^d^'^^'W^^^' 

which determines -j-j : 

da 
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da* da da' da * 

which determines , ,^ : 

oaaa 

da'* dada' da" 

which determines ,-L* ^ and thus all the dififerential coeffi- 

da 

cients of G. may be determined. 

In like manner all the succeasive differential coefficients 
of JD may be calculated. 

We are now in a position to simplify the expression for 
F. We have (Art. 47) 

^4V da' da'*/ 

■^4-(^^'-^^'^f-^2'^''^'+««'^)-'''°^(— ') 

+ 7 aa'i>j tan i tan i' cos (ft — fl') + ... h . 

The following proposition will be found useful. 

d*C d*C 

58. To shew that 3— r%= — D,, and that ^ — rS-= - D*. 

da da * da da * 

We have 

- G^+ C^cos(f>+G^coH2<l> + ... = {a+a'^-2aa'co3^y^ 

^, . \dG. dG, . , dO, ^ . , 

therefore ^ -j-^ + j- cos 6 + -r-^ cos 20 + . . . 

z da da ^ da ^ 

= — {a'-a' cos (f>) (a* + a'* — 2aa' cos ^)""' ; 



-I 
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therefore 

4- 3 (a — a' cos ^) (a' — a cos <^) (a* + a" - 2aa' cos <^)"^ 
= cos <f> (a* + a'* — 2aa' cos <^)"^ 
+ 3 [aa (1 + cos' <f>) - (a' + a'*) cos <^} (a' + a ' - 2aa' cos ^) 
= cos <f> {a^ + a'' — 2aa cos ^)"^ 
+ 3 [ad' sin'^ ^ — cos ^ (a* + a'^ — 2aa' cos <^)} 

(a* + a ' — 2(ya' cos ^)"' 
= — 2 cos ^ (a* + a'^ — 2aa' cos <^)"^ 

+ Saa' sin' ^ (a* -{- a" - 2aa cos ^)"^. 

8 1 

Now (a' + a'* - 2aa' cos ^)"^ = ^ Z)^ + A ^^^ <A + • • • ; 

diflferentiating with respect to 

3aa' sin <^ (a' -fa'"— 2aa' cos <^)~5 

= Dj sin <^ 4- 2i>3 sin 2<^ + . . . ; 

therefore ^ ;j-±, + ;t-^, cos-^ + . . . 



2 cfa oSa (2a da' 

= — 2 cos 



^fg A + A^^s^+ ...j 



+ sin ^ (2)j sin ^ + 22), sin 2^ -f- ...)> 

whence, equating the parts independent of <f>, and also the 
coefficients of cos <f> 

1 (fg. _ Z),_ J. 

2rfado'~ ' 2 ~ 2 ' 

tfq n. 
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dC 

59. Since —r^ is a homogeneous function of a and a' of 

- 2 dimensions, 

a— » + a'-^^=-2— »• 
da^ da da da ' 

,, - dG, a'd'C, aa' d^G, 

therefore o.-j-^-k--^ -j-r = - -«- j— tt 

da 2 aa 2 da da 

= ^ aa' Dj. 

Similarly. a'§ + ^''^" = |aa'A. 

Hence the coeflScients of e" and e'* in the expression for F 
are each equal to ^ oa' D . 

o 

Again, since C^ is a homogeneous function of a and a' of 

— 1 dimensions, 

"d^^"!^' — ^^' 
hence the coefficient of ee' cos (to- — -or') 

= i(2C,-aa'D.); 

but (Art. 55) 2AC; = aa' {D^^ - D^,); 

therefore, making k = l, 

2C. = aa'(Z),-2>,); 
hence the coeflScient of ee' cos {w — w') 



= — -T aa B,. 
4 



Again, (Art. 66) 

dC. 



^^-aD. + a^D,, 
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59 



theiefore 






Similarly, d'D, + a' ^^ = axiD, 



da' 

Hence the coefficients of tan* i and tan' i' are each equal 
to — - aa'D., 

o 

r 
»• V 

Fiaally, the expression for F becomes 
i^= to' I ^-^ + 1 oa' 2>^ (e* + e-^ - i aa' D, ce' cos (isr - w') 

— 5 aa'D^ (tan' i + tan* i) + ^ aa'D^ tan i tan i' cos (Q - ft') 



8 



4 



+ ...K 



CHAPTER IV. 

SECULAR VARIATIONS OF THE ELEMENTS OF THE ORBIT. 
STABILITY OF THE PLANETARY SYSTEM. 

60. We have seen in the preceding Chapter, that the 
disturbing function, when developed, consists of two parts; 
the one independent of the time explicitly, the other in- 
volving it under a periodical form : we shall consider sepa- 
rately the eflfects of these two parts. In the present Chapter 
our attention will be directed to the first or non-periodical 
part of iJ, which we have denoted by F, The inequalities 
thus produced in the elements of the orbit are termed «6cw?ar, 
in consequence of their very slow variation. 

61. By differentiating the expression for F in Art. 59, 
with respect to the elements, we obtain 

de "' 

dF m! ,r\ , • / f\ 

=:— aa l)^ee sm («r — -cr ), 



d 



nr 



dF m f T. m ,j^ , . , 

^=:—.aa V^e -- -r-CLa B^e cos («r — "or ), 

■=^ = — -Y^ aa^D^ tan i tan i sin (fl — fl'). 
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^ = - ^- aa jDj tan z + ^ aaD^ tan % cos (ft - fl ), 

-3— = an expression similar to F, 

Substituting these in the formulae of Art. 39, and neg- 
lecting powers of e, e\ tan i, tan i higher than the second, 
we have 

dt ' 

de m'mnJ'a' r, , . , ^ 

di^ 4^ A^ sin (t^ - t^rO, 

dt ^ 4ifM A^ani sm (ft - il ), 



tan I 



T — i^i {tan z*— tan i' cos (ft — ft')}, 



+ A^ee' cos (-bj — -or') + J.^ tan i tan { cos (ft — ft'), 

where in the last expression, -4, -4^, &c., have been written 
to denote certain functions of a and a\ 

62. To calculate approodmately the secular variations of 
the elements of a pUmet's orbit, in a given time. 

Let a^, e^y vx^, &c., be the values of the elements at some 
given epoch; a^ + Ba*, e^-h Be, w^-hBrir, &c., their values after 
an interval t: then Ba, Be, Bw, &c., are the required variations. 
By Maclaurin's Theorem, 

* It wiU be shewn in Art. 64 that Za is always zero« 
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^■■(w).-($)|-- 

• • • ^~" • • • 

which may be carried to any required degree of accuracy, but 
in practice the first two terms will generally be sufficient. 

We have supposed the variations of the elements required 
at a time t after the epoch ; if they be required at a time t 
before the epoch, we have only to change the sign of t in the 
above. 

We may remark that (-7-] , (-^^r) , &c. are of the order 
of the disturbing force, since they involve the first power of 
^^' • [j^] * {~^) ' ^^' ^^^ oftli^ second order; for, since the 

expressions for jt, -^^ , &c., involve elements, their differ- 
ential coefficients will involve the differential coefficients of 
those elements, and thus, by substitution, W will be intro- 
duced: similarly (-tjz) > [~j^) > *^c., are of the third order, 

and so on. 

In the short period of one year all terms after the first 
may be neglected, so that putting ^ = 1, we have 

& = (I); &c. 

Hence the coefficient of t in the above formulae is called 
the annual variation, 

63. Since the elements of the planetary orbits are con- 
tinually changing, it will be interesting to shew that the 
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dimensions of these orbits, and their inclinations to the 
ecliptic, nevertheless fluctuate between very narrow limits. 
This constitutes what is termed the Stability of the Plane- 
tary System : in order to establish it, it will be necessary 
to prove the stability (i) of the mean distances, (ii) of the 
excentricities, (iii) of the inclinations. 

64. To prove the stability of the mean distances of the 
planets from the Sun, and of their mean motions. 

da 
By Art. 61 -^r = 0, so that a is constant. Now it will be 

shewn in a subsequent Chapter (see Art. 91), that to the first 
order of the disturbing force, the periodical terms of R can 
produce only periodical variations; consequently, to this 
order, the mean distance is susceptible of no permanent 
change*. The same is true of the mean motion n, since 

J 
it = — j , and fju does not alter. We are hereby assured of the 

a 
impossibility of any of the bodies of our system ever leaving it 
in consequence of the disturbances it may experience from 
the other bodies ; and this secures the general permanence of 
the whole, by keeping the mean distances and periodic times 
perpetually fluctuating between certain limits (very restricted 
ones) which they can never exceed or fall short of. 

This result may easily be extended to all orders of the 
excentricities and inclinations : for since nt-\- € always occurs 
in R as one symbol, e cannot occur in F because t does not, 

so that — , and therefore -j- is zero. 
lie at 

* This result is also tine when the square of the disturbing force is 
included: for the demonstration the reader is referred to Pont6coulant's 
Sysi^me du Monde, Tome I. p. 395 (2nd edit.). 
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65. To prove the stability of the excentricities of the pla- 
netary orbits. 

We will first consider the case of two planets only. By 

Art. 61, 

de m'na^a* T^ / • / /\ 

^. ., , de* mna'^a -r,, * f » v 
Sinnlarly, -77- = ^r -^a ^ sin («r — -or). 

Now since 2), is the coefficient of cos 2<f> in the develop- 
ment of 

(a' + a * - 2aa' cos ^)"*, 

an expression in which a and a' are similarly involved, it 
follows that 

Hence, multiplying the above equations by — e, -7-, e', 

respectively, and adding, we have 

m de , wl ,de ,. 
na at na at 

therefore, since a experiences no secular variation, 

na na 
A similar equation holds for any number of planets. Re- 
placing for convenience -7 — - by (a, a'), we have 

— = — mna (a, a') e' sin («r — -bj') 

— m"wa (a, a'') e" sin (tar — -bj") — . . . 

^-= — TnnV (a', a) 6 sin (w' — -or) 

-. mVa' (a', a") e'' sin (tsr' - «r") - . . . 
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dt 



>• • •. 



= — mti'd' {d\ a) t sin (ot"— w) 



— mna (ajajesin^-cr — -or) — .. 



Since D^^D^, it follows that (a, a') = (a', a): hence mul- 
tiplying these equations by — e, -7-; e\ &c., and adding, 
we obtain 



m de mf , de' , w!' ,, de" ^ _^^ 



whence by integration, 

2 ( — e* ) = const. 
\na / 

Since wa = /- , this may be written 

S ( w V^*) = C^p 

Now observation shews that all the planets revolve round 
the Sun in the same direction, so that the mean motions w, n\ 
vl\ &c. are of uniform sign. Hence all the terms of the left- 
hand member of the above equation are positive. 

We learn also from observation that the excentricities of 
the planetary orbits are at present very small indeed, and 
in the case of the Asteroids, the masses are very small^ 
Hence the constant must be small. Since, then, all the terms 
of the first side of the equation are positive, and their sum 
always equals a small constant, it follows that every term is 
small ; and therefore, except in the case of planets of small 
mass, such as Mercury, Mars, J uno, &c., that the excentricities 
must remain permanently small. 

The stability of the excentricities, however, is not con- 
fined to the larger planets ; we shall give another proof of 
c. P, T. 5 
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this important theorem and that of Art. 66 in the following 
Chapter. 

66. To prove the stability of the inclinatioTj^ of the planes 
of the planetary orbits. 

By Art. 61 -r, = — r D^ tan i sin (il — II'). 

Similarly, -^ =? — j-^ — 2)/ tan i sin- (il' — II). 

As in the preceding Article, it may be shewn that 
D^' = Dj. Hence, multiplying the above equations by 

m . . rn! , „ 
-^— tan %y —J—, tan % , 
na na 

respectively, and adding, we have 

w ^ ,di , m . ,,di r, 
— tant-T- + -7-7tan* -:t: = 0, 
na at • na at 

or to the same order of approximation, 

m ^ .d(tant) m' ., d(tani') ^ 

— tani— — — -4 -7— , tani ^ ,, — ^=0; 
na at na at 

therefore — tan' t + -7-7 tan* i' == const., 

na na 

or, since na= a/- , 

m hja tan'i + m tsja' tanV =* C, 

A similar equation would (as in the case of the excentri- 
cities) be true for any number of planets. Now the inclina- 
tions of the planetary orbits to the ecliptic are at present very 
small; hence, if we take for our fixed plane of reference a 
plane coinciding with the present position of the ecliptic, and 
except the case of planets of dmall mass, it follows, as in 
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Art. 65, that their inclinations to thia plane must always 
remain very small. 

67. It should be noticed that the assumption in the two 
preceding Articles, that ji is the same for all the planets, is 
equivalent to neglecting the square of the disturbing force : 
for, let fi refer to the planet m and jj! to m'; then, if M 
denote the Sun's mass, ^=^M+m, ij!^M'\-m\ so that ii 
differs from /* by a quantity of the order of the disturbing 
force : since, then, the expressions in which fx and fi' occur 
are themselves of the first order, it follows that the error 
introduced by supposing jx and fi equal is of the second 
order. 

68. The results of Arts. 65 and 66 may also be obtained 
directly from the equations of motion. We have (Arts. 20 
and 15) 

dt dd de dvr ' 

or, replacing R by F, since we have seen that -r- = 0, 

dh^dF 
dt dxsr 

Multiplying both sides of this equation by m, forming 
similar equations for each planet of the system, and adding, 
we have 

^(»S)-^(»S)- 

Now on referring to the expression for F in Art. 59, since 
VT occurs only in the term " — da' D^ee cos {^ --vr'), it is 

(dW\ 
m^— j = Oj hence our equation becomes 

S(»g).0, 

5—2 
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whence by integration 

2(wA) = const.; 

or since Ji? — /jLa(l— e*) = (Jlf + m) a (1 — e^, 

% \mt^(Ma) (l +-^^ + ...]\1 -^ -{- ...U = const. 

Since a is constant as regards secular variations, and our 
approximation extends only to the second order of the excen- 
tricities, we have to the first order of the disturbing force 

S {m *Jae^) = C, 

the equation of Art 65. 

Again, by referring the motion to the fixed plane of refer- 
ence, we obtain 

d .J ^ dR dR dR dR .. -^v 
^^(Acost) = ^ = ^- + ^- + ^(Axt.l6), 

and considering the whole planetary system, we have as 
before 



23 






Now on referring to the expression for jP, it is easily seen 
from the forms imder which to- and H occur, that 

ft 

hence our equation becomes 



2 jm ^ (A cos in = 0, 



whence by integration 

2 {mh cos i) = const*, 

* It should be borne in mind that this and the equation Z (mA)= C, which 
are those which would be obtained bj conservation of areas were it lawful to 
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or S j mA f 1 — 2 sin'' k ) f = const., 

which, from above, may be written 

2 (mh sin* ^ j = const.> 

and proceeding as before, we obtain to the second order of 
the excentricities and inclinations, and the first of the dis- 
turbing force 

2 ( m tja sin* ^ j = const., 

or, since to this order of approximation, the inclinations, their 
sines and tangents are equal, this may be written 

2 (m ^a tan' i) = C, 

the equation of Art. 66. 

We may remark that the equation 

2 {nth cos i) = const. 

is of itself sufficient to establish the stability both of the 
excentricities and inclinations. For, proceeding as before, 
we obtain 

s{mV(ifa)(l+|j+...)(l-^+...)(l-f + ...)}=con3t., 

which, to the second order of the excentricities and inclina- 
tions and the first of the disturbing force, gives 

2 {m hjai) + 2 (m V^t**) = ^> 
or to the same order of approximation 

2 (m Vae') + 2 (w V^t tan* i) = C 

Since we know from observation that all the planets 
revolve round the Sun in the same direction, all the radicals 

assume the principle, are tme only as regards secular variations. To 
assume their actual truth would be to neglect terms in B. due to the Sun's 
motion, of the first order of the disturbing force. 
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in this equation must be taken with the same sign. Also, 
since the excentricities and inclinations are at present very 
small, the constant must be small. Hence it follows, as in 
Arts. 65 and 66, that the excentricities and inclinations must 
always remain very small. 

69. From the preceding analysis we draw the following 
remarkable conclusion : The fact that the planets revolve about 
the Sun in the same direction, ensures the stability of the 
planetary system. The converse of this would not necessa- 
rily be true, as we shall see in Art. 75 : the numerical rela- 
tions of the dimensions and positions of the orbits of the 
planets, might be such as to ensure stability, although they 
revolved in opposite directions. But the above is indepen- 
dent of particular numerical relations. 

70. The results at which we have arrived with regard to 
the stability of the planetary system are of especial interest. 
In consequence of the changes in the elements it might have 
been supposed that the orbits would ultimately undergo 
such alterations in their dimensions as to bring the planets 
into collision or hurry them into boundless space. Or even if 
no such violent catastrophe occurred, a derangement of the 
seasons might seriously have interfered with the physical 
comfort of man*. But our analysis shews, (and the results 
are confirmed when the approximation is carried further,) 
that in so far as the mutual attractions of the Sun and 
planets are concerned, the dimensions and position of the 
orbits will for ages remain nearly the same as they are at 
present, i.e. nearly circular in form, and but little inclined to 
each other, thus aflfording a beautiful illustration of Gen. viii. 
22: "While the earth remaineth, seed-time and harvest, 
and cold and heat, and summer and winter, and day and 
night shall not cease." 

* Se^ Herechers Outlines of Astronomy* 



CHAPTER V. 

r 

SECULAR VARIATIONS OF THE ELEMENTS CONTINUED. 
INTEGRATION OF THE DIFFERENTIAL EQUATIONS. 

n 

71. In Art. 62 we have given a method of calculating 
the secular variations sufficiently accurate for the practical 
purposes of astronomy, but in order to understand their real 
nature, and thus to examine more fully into the important 
question of the stability of the excentricities and inclina- 
tions, it is necessary to proceed to the actual integration of 
the equations of Art. 61. This we are enabled to do by a 
method due to Lagrange. 

72. To integrate the equations for the eccentricity and 
longitude of perihelion. 

We have (Art. 61) for the planet m 
with similar equations for the planet m\ 
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We shall be able to reduce these to a system of linear 
differential equations if we assume 

ti = ^ sin cr, V ^e cos«r, 

u'= c'sin m, v' = e' cos «•' ; 

therefore 3- = 6 cos ta- -j~ + sm ta- 3- » 

at at at 

Substituting the values of -1- and -y- j and writing a 



dt dt 



J, Tnfna^a , 

for — 2 > we nave 

4/i 



mnfa'^a 



3- = a (-Oj^ cos «r — JDje' cos'cr') 

Similarly, ^ ~ ^ (A^' "" A^)* 

In like manner for the planet m', writing a for , 

we have 

|^ = a'(A«-A«'). 

The forms of these equations suggest the following par- 
ticular integrals : 

u = M sin (fft + y), v =M coB(gt + y), 
u'=M' sin (gt + 7), v' = M' cos {gt + 7). 

Substituting these in the differential equations, we obtain 
from either of the first two 

gM^a{D^M^D^M% 

and from either of the last two 
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eliminating the ratio M \ M 

or /-(a + a')A5' + «a'(A'-A') = 0: 

and the roots of this equation will be real and unequal, real 
and equal, or impossible, according as 

(a + aj D^ - 4aa' {D^ - D^) 
is positive, zero, or negative. Now 

(a + aTD,^ - 4aa'(A' " ^D = (« " «? A' + 4^a A*> 
A positive quantity, since w, n' and therefore a, a are of like 
sign. Hence the values of ^r will be real and unequal: denote 
them by g^,g^, and let 7^,7^; Jf,, M^\ if/, M^\ be the cor- 
responding values of 7, My M' respectively. Then the com- 
plete solution of the differential equations will be 

w = M^ sin {gf + 7,) + if, sin {gj; + 7,), 
v = if, cos {g^t + 7j) + M^ cos (gr^i + 7,), 
u' =M; sin (gr^^ + 7^) + if/ sin {gji + 7,), 
t;' = if; cos {g^t + 7j + M; cos (^r^^H- 7,). 

Of the constants in these equations, four are arbitrary and 
must be determined from observation. We have 

e« = ti« + t;« = if,^ + if,' + 2if,if,cos{(5r,-<7,)^ + 7i-7j* 

V M^ cos {g^t + 7,) + if, cos {gji + 7,) ' 
with similar equations for e' and -cr'. 

73. Had we considered a system of several planets, we 
should have obtained by a like process 

+ 2M^M^ cos ((^r,- 5- j < + 7i- 7,1 + • • • 
^^^ ^ Jf, sin (.y,<+7,) + Ji;8iD (gr,< 4- 7,)+.af,8in(.y,<+ 7,)+ . . ■ 

if, cos (^,<+7,) +i/,COS (gr^e +7^ + J/,COS(^,<+ 7,)+ . . . ' 

with similar equations for each of the other planets. 
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74. We may faence infer the stability of the excen- 
tricitifes. From the form of the expression for e it appears 
that e' cannot be greater than 

and therefore that e cannot exceed 

-Jf^ + ifj + JIfjH-..., 
these quantities being all taken with the same sign. 

Thus, by determining the numerical values of Jf^, M^^ 
&c., for any particular planet, we may assign an actual limit 
which its excentricity can never exceed. For the principal 
planets, Jf^, M^, &c., are found to be exceedingly small, so 
that the excentricity must always remain veiy small. 

If we confine our attention to two planets, we see 
from the expression for e in Art. 72, that the excentricity 
fluctuates between the limits M^ + M^ and M^-^M^, The 

period of these chancres = , and is the same for each 

planet : it appears also from the equation 

obtained in Art. 65, that the maximum of each excentricity 
takes place at the time of the minimum of the other. 

As an illustration, take the case of Jupiter and Saturn. 
Sir John Herschel finds that 

^, = 21"-9905, ^,=:3"-5851; 

Jf, = - -01715, M^ = -04321, for Jupiter ; 

Jf/= -04877, Jf/= -03532, for Saturn: 

the year 1700 being taken as the epoch*. Thus we obtain 
for the greatest and least excentricities that Jupiter's orbit 
can attain, -06036 and '02606 respectively, and for those of 

* Article Physical Astronomy in the Encyclopadia Metropolitana. 



SECULAR VARIATIONS. 75 

Saturn, •08409 and •01345, quantities exceedingly minute; 
while the period of these changes amounts to no less than 
70414 years. 

75. It appears from the preceding Article that the sta^ 
bility of the excentricities is a consequence of the periodical 
form of the solution of the differential equations, a result 
which depends upon the fact that g^ and g^ are real and 
unequal. Now we have seen that in order that this may be 
the case, it is only necessary that 

(a+a7i),^-4aa'(i),*-Z),0 

shall be positive, a condition which might be satisfied if the 
signs of 71, Thy and therefore of a, a! were different. In this 
case, then, the stability would still subsist. Let us, however, 
consider what would be the effect of equal or impossible roots 
to the quadratic from which g is found. In the former case 
a term would be introduced into w, u\ v, and v' proportional 
to the time, and in the latter the periodical terms would be 
replaced by exponentials. Consequently the excentricities 
would increase indefinitely with the time, and the stability 
would no longer subsist. 

76. We now proceed to examine the expression which 
has been obtained in Art. 72, for the longitude of perihe- 
lion, viz. 

M^ cos {g^t+ 7,)+if, cos(5r,^ +7,) ' 

. d^ ^ gA'+g2^,'+ (9i+ff,) M^M.cos {(g.-g,) t+ 7^- 7«} 
•• dt M^'+M,' + 2M^M,cos[{g,-g,)t + y,--y,} ' 

The maxima and minima values of ta-, if such exist, will be 
found by equating -jr to zero. Thus 



dt 



8 
3 
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If this (disregarding sign) be not greater than unity, the 
perihelion will oscillate, the period of a complete oscillation 

beinof the same as that of the excentricities, viz. : but 

if, as is the case with Jupiter and Saturn, this be greater 
than unity, the longitude of perihelion has no maximum or 
minimum, and the perihelion moves constantly in one direc- 
tion. 

Again, 

d^ 1 fa-fir,)(Jf.'-3f,') 1. . 

dt 2 M*^M:+ m,M^ COS {{g,-g;) t + 7^- 7,} ^ 2 ^^' ^, ^^^ 

Hence when e is a maximum or minimum, -rr will be 

dt 

either a maximum or minimum, and the apsidal line will be 
moving most rapidly or most slowly, diflferent cases occur- 
ring according to the signs and magnitudes of the quantities 
involved. 



77. Wheal the apsidal line oscillates^ to find the extent 
and periods of its oscillations* 

We have (Art. 72) 

^^^^ -^iSin(gr,^-h7,) + 3faSin(ff,^ + 7,) . 
M^ cos {g^t + 7,) -h if, cos {gj, + 7,) ' 

., ^ . , ^ V tan'Br-tan(5ri + 7) 

therefore tan (w — g.t — 7.) = z — 7 7 — ^f-S — ^ 

\ if I /I/ 1 4. tan -cr tan (^r^t -H 7i) 

_ — -JfjSin-^ 

"" Jlf J + J/j, cos >/r ' 
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Also by the last Article, if r be the least positive angle 

, . . g.M*-\-gJd^ 

whose cosine is p — * — ririiVi 

(5^1 + 5^2) ^A 

e^^= (fft + g^ M^M^ (cos f + cos t), 

Diflferent cases will occur according to the signs of if,, 
Jfj, &c. Suppose M^y M^ of like sign, g^ and g^ positive, and 

g^ greater than g^. Then -^ increases as t increases, and -77 

will be negative, or the apsidal line will regrede, while 
cos -^4- COST is negative, i.e. so long as -^ is between 

(2/1 — 1) TT — T and (2n — 1) tt + t; -^- will be positive, or the 

apsidal line will progrede, while -^ is between {2n — 1) tt + t 
and (2n + 1) tt — r. 

To find the angle through which the apsidal line regredes 
and the period of the regression. Let t\ f be the values 
of f, «r', tsr" the values of -cr corresponding to the values 
(2/1 — 1) TT — T and (2/i — 1) tt + t of -^ : then 

(5'i-fl'2)«" + 7i-7, = (2n-l)7r + T, 

. , , .f V —if. sin T 
tan ('cr — at t — ry ) = -— J- 

. / // y/ \ if. sin T 

tan(OT — flrf — 7.) = iri^ — S^F , 

\ H\ 11/ if^-if^cOST 

From these equations the values of t\ t'\ ts\ «•", may be 
found, and thus «■' — -bj'' the amount of regression will be 
known. The period of regression 



78 
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In like manner the amount and period of the progression 

2 (tt — t) 
may be obtained. The latter will be found to be — ^-- -. 

The period of a complete oscillation will be the sum of 
the periods of the regression and progression, that is ^^ 



which agrees with the preceding Article. 



9i-9i 



78. The motion of the centre of the instantaneous ellipse 
in consequence of the secular variations of e and vr may be 
exhibited geometrically as follows. 

We have, by Art. 72, 

6 cos lar = M^ cos {gji + 7^) + M^ cos {gji 4- 7a), 
e sin -cr = M^ sin {g^ + 7^) + M^ sin {gj, + 7^^. 

Let a circle be described in the plane of the orbit with 
its centre 8 coinciding with that of the Sun, and its radius 
equal to M^a^ where a is the mean distance. Let a point P 
describe this circle uniformly with a velocity g^^ starting from 
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0. Agaiq, with centre P and radius equal to M^d let a^othe^ 
circle be described, and let a point Q describe this circle 
uniformly with a velocity g^, starting from G, Let SL be 
the line from which longitudes are reckoned, and draw PK 
parallel to it: then if the angle OfiW be equal to 7,, and 
GPK to 7j, the angle P8N will be equal to gf + 7^, and 
QPK to 5^2^ + 7j. Produce QP to meet the circle again in i2, 
and draw QN perpendicular to SN. Then, supposing M^^ and 
M^ to be both positive, we have 

SN = SP cos PSN + PQ cos QPK 

= M^a cos {g^t + 7,) + M^a cos {g^t + 7,) 
= ae cos -cr. 

Similarly, it may be shewn that 

QiV=aesin'cr. 

Hence, the apse being supposed to move from L in the 
direction contrary to that of the hands of a watch, Q will be 
the centre of the instantaneous ellipse. 

If J/j be positive and M^ negative, it may be shewn in 
like manner that the centre of the ellipse will be R. If 
Ifj, il/j be both negative, join Q>S and produce it to Q so 
that SQ' = SQ : then the centre of the ellipse will be Q', 

A similar construction will of course apply for the motion 
of the further focus. 

It is easily seen from the above that the excentricity is 
least when Q is in the line SP and greatest when Q is in the 
line SP produced. Hence the maximum and minimum 
values of the excentricity are ifj + Jfef, and M^-^M^ respec- 
tively. (Art. 74.) 
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79, To integrate the eqiuitions for the inclination and lanr 
gitude of the node. 

We have (Art. 61) for the planet m 
di mncfa' 



dt 4ifA 



D^ tan i' sin (fl - fl'), 



tan * -^ = T -Pj {tan % — tan i cos (XI — ll )} ; 

with similar equations for the planet tw'. 

To integrate these, assume 

p — tan i sin fl, 5^ = tan i cos fl, 
p' = tan % sin 11', g' = tan H cos XI' ; 

therefore -^ = tan i cos XI -rr + sin XI (1.+ tan'i) -7: , 
dt dt at 

Substituting the expressions for -^ and -^ , and writing 

a for — 7 , since tan^z -^ being of the third order may be 

omitted, we have 

-^ = olD^ (tan i' cos XI' — tan i cos XI) 
Similarly, ^ = a^i(p-/)- 



Also for the planet m', writing a' for 



mn'a^a 



4/i ' 
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The forms of these equations suggest the foUowing par- 
ticular integrals : 

jj = JVsin (ht + S), q = JVcos (ht + S), 

pr^N'sin Qvt + S), ^^N'cos (ht+ S). 

Substituting these in the differential equations, we obtain 
from either of the first two, 

and from either of the last two, 

eliminating the ratio N : N\ 

(h + ojDJ (h + a'D,) = aa'i),«, 
or h*-\'(a + a')DJi = 0; 

therefore A = — (a + a') D^, or A = 0. 

Denote the former by h^, and let 8„ S^, JV,, N^, J\r', iV/, 
be the values of S, N, N' corresponding to A = A^ and A = 0. 
Then N^ =N^, and the complete solution of the differential 
equations will be 

p = iV; sin (h^t + S^) + iVT sin S^, 

gr = i\r, cos (A,« + 8,) + N^ cos S,, 

/= JV; sin (A,* + S,) + N, sin S,, 

2' = JV/ cos (A,* + SJ + JV, cos 8,. 

Of the constants in these equations, four are arbitrary, 
and must be determined from the known yalues of i and fl at 
some given epoch. 

We have then 

tan« i =/ + 3« = N,' + N,' + 2N'^N^ cos (hj; + 8, - S,\ 

t^^ fl ^i> -, ^t sin (A,^ + 8,) + JT, sin 8, 
q JN'^ cos (A^^ + Sj) + i^2 cos 8^ * 

with similar equations for i' and D.\ 

C. P. T. 6 
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Had we considered a system of several planets, we should 
have obtained a result precisely similar to that of Art. 73. 

80. From the form of the expression for tan't, the sta- 
bility of the inclinations may be inferred. For it may be 
shewn, as in Art. 74, that tan % can never exceed 

these quantities being taken with the same sign. Since 
then N^, JV^, N^ &c. are found to be very small, it follows 
that the inclinations must always remain exceedingly small. 

In the case of two mutually disturbing planets, we learn 
from the expression in Art. 79, that tan i fluctuates between: 
the limits iV^^ + iVj and N^-^ N^. The periods of these 

changes are the same for the two planets, being —-7- ; and 

as appears from the equation of Art. 66, the maximum of 
each inclination will take place at the time of the minimum 
of the other. 

In the case of Jupiter and Saturn, the period is 50673 
years ; the maximum and minimum inclinations of Jupiter's 
orbit to the ecliptic are 2* 2' 30" and 1" 17' 10", those of 
Saturn's orbit 2" 32' 40" and 0° 47'. 

81. "We now proceed to examine the expression which 
has been obtained in Art. 79 for the longitude of the node. 
We have 

.^^ Q _ -^t sin (h,t + B,) + N, sin 8, . 
^^ ~ 2f, cos (h^t + S,) + N^ cos 8, • 

The maxima and minima values of £1, if such exist, will 
be found by equating -t- to zero. Thus 
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If this (disregarding sign) be not greater than unity, the 

node will oscillate, the period of a complete oscillation being 

27r 
the same as that of the inclinations, viz. —-7- . But if it be 

greater than unity, there cannot be any stationary positions, 
and the node will move continually in one direction. 

It may be shewn, as in Art. 76, that the motion of the 
node will be fastest or slowest whenever the inclination is 
either a maximum or a minimum. 

82. When the line of nodes oscillates, to find the extent 
and periods of its oscillations. 

It may be shewn as in Art. 77, that if ^jt be written for 

Aji+8j — Sjj, and t denote the least positive angle whose 

. . K 
cosme IS -^^ , 



tan (n - J = T.J /,^ — ^—7 

N^ + iV^ COS Y 

cos T sin yjt 



1 + cos T cos yjt ' 
and tan* * -ir = h^N'^N^ (cos yjr + cos t). 

Different cases will occur according to the signs of iV^, 
JVj and Aj. Suppose iV^, N^ of like sign, and h^ positive: 
then ylr increases as t increases, and the line of nodes re- 
gredes so long as -^ is between (2n— 1) tt — t and {2n — 1) 7r+ t, 
and progredes so long as '^ is between (2n — 1) tt + t and 
(2w+l)7r-T. 

Let fl', 11" be the values of fl corresponding to the values 
(2n — 1) TT - T and (2n — 1) tt + t of -^ ; then 

tan (fl' — Sj) = cot T, 

6—2 



qgH)p«P^^(9a«pi>va«spav^v^Hpm«willppi> ■ ni I ■ 1 nri V— ^fv'" •. »"^»r"""T^"W9^ 
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therefore fi' — S^ = m^r + ^ — t, 

n"-S^ = (m-l)7r + | + T; 

therefore fl'- X2" = 7r-2T, 

which is the angle through which the line of nodes regredes. 
Also the period of this regression may be shewn as in Art. 77 

2t 

to be T-. Similarly, the angle through which the line of 

nodes progredes may be shewn to be tt — 2t, and the period 

i^ *u • 2 (tt — t) 

of the progression -^-r — - . 

The period of a complete oscillation will be the sum of 
the periods of the regression and progression, that is -r- , 
which agrees with the preceding Article. 

The remaining cases corresponding to diflferent arrange- 
ments of the signs of iV^, N^ and \ may be treated in like 
manner. 

The mean value of 12 is lYiir -f 8^, whatever be the signs 
of N^y N^ and \, and the mean value coincides with the 
true whenever sin -^ = 0. Since, then, '^ is the same both 
for the disturbed and disturbing planet, the nodes of both 
orbits will arrive simultaneously at their mean positions. 

In the case of Jupiter and Saturn, N^ is for each planet 
numerically less than iV^, so that the node oscillates ; the 
extent of oscillation being 13** 9' 40" in Jupiter's orbit, and 
31** 56' 20" in that of Saturn on either side of their mean 
position, the ecliptic being taken for the plane of reference, 
and supposed immoveable. 
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83. To shew that the inclination of the orbits of two 
mutually disturbing planets to each other is approximately* 
constant. 

If 7 denote this inclination, we have by Spherical Trigo-* 
nometry, 

cos 7 = cos i cos i + sin i sin i' cos (11 — fl') 

= cos i cos t' {1 + tan i tan % cos (fl — Si')] 

= (1 + tan't)-i (1 + tan'*')-* {1 + tan i tan % cos (ft - ft^} 

1 
= 1 — ^ {tan't + tan'i' — 2 tan { tan % cos (ft — ft')}, 

if we neglect small quantities of orders higher than the 
second. 

Now tan'i + tan'i' — 2 tan i tan i' cos (ft — ft') 
= p'+3'+i)'' + 2'"-2(|)p' + j2') 

therefore 1 — cos 7 = ^ {N^^ — N^f, 

whence it follows that 7 is constant. 

84. The equations which give the secular variations 
of the node and inclination may be explained geometrically 
as follows*. 

The equations to be interpreted are 

^ = i^, sin (A,< + 8,) + i^, sin S„ 

g = -N'j cos {h^t + 8j) + N^ cos Sj, 

where p = tanisinft, j = tanico8ft. 

^ This elegant geametrieal explanation is due to Mr H. M. Taylor, M.A., 
Fellow and Tutor of Trinity College, Cambridge. Oxford^ Cambridge and 
DMin Messenger of Mathematics, Vol. m. p. 189. 
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Let NA^ NB be the intersections of the planes of refer- 
ence and of the orbit respectively with a sphere of radius 
unity, the centre of which coincides with that of the Sun : 
ABZP the great circle of which the pole is N: Zand Pthe 
poles of the great circles NA^ NB : L the point from which 
longitude is measured 

Then LN^ fl, and PZ= BA^i 




Now project PZ and the other great circles by radii 
drawn from the centre of the sphere on the tangent plane at 
Z: then tantsinfl and tantcosfl are the Cartesian co- 
ordinates of the projection of P referred to the projection 
of LZ as axis of y^ and a line at right angles to it as axis 
of cc. 

Suppose P' the projection of P, then if the co-ordinates 
of P' be X and y, we have 

x^N^ cos (hf + Sj) + N^ cos Sj, 
y = Jf, sin {\t + S,) + N^ sin S,. 

These equations shew us that P always lies on a circle of 
which the centre is at the fixed point {N^cosi^^ j^T^sinSj), 
and the radius is N^: also that P' describes this circle with a 
uniform angular velocity h^. 
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87 




We may hence arrive geometrically at the results of 
Arts. 81 and 82. If C be the centre of the circle, we have 

angle CZx = 8, , angle PVO = h^t + S, - S, ; 

and if P'Z be joined, 

P'Z=i, angleP'Za? = a 

Suppose JV^j, N^ both positive ; and first, let N^ be less 
than N^, 

Now the only time when the node will be stationary will 
be when P' is moving directly towards or directly from Z, 
that is at jT, T' the points of contact of tangents to the circle 
drawn from Z. As P' moves from T to T, fl increases, or 
the node progredes, and as P" moves from Tto T\ 12 decreases, 
or the node regredes. 

If, then, T denote the angle TCZ I cos"^ W 1 , the period of 

the progression will be the time P' takes to move from T' to 

T, that is —^ ; the period of the regression will be the 

2t 
time P takes to move from T to T, that is -r-; and the node 
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is stationary whenever \t + 8^ — 8^ = (2n + 1) tt ± t. Also 
the angle through which the node oscillates is TZT' which 

= TT - 2t. 

Secondly, suppose N^^N^, Then Z will lie on the cir- 
cumference of the circle which P' traces out In this case T 
and T coincide, and the no^e after becoming stationary 
begins to move in the same direction as before. 

Thirdly, suppose N^ greater than N^. Then Z lies within 
the circle, and no tangents can be drawn from it to the circle: 
from this we see that the node is never stationary^ 

It is easily seen that in all three cases the maximum 
and minimum values of the inclination are ZO and Zffy that 
is -ftTj+^^and N^-^ N^ respectively. And in all cases the 
node moves fastest or slowest when P' coincides with or 0', 
that is, whenever the inclination is a maximum or minimum. 

The above geometrical construction also affords a proof of 
the theorem of Art. 83. 

Since in the case of two mutually disturbing planets the 
quantities \, 8^, N^j S, are the same for both, it follows that 
the point P'' for the second planet traces out a circle con- 
centric with that traced out by P', and with the same uni- 
form angular velocity A^; also that the two points P*, P" 
always lie in a common radius vector through G. Now the 
angle subtended by P'P' at the centre of the sphere is the 
inclination of the two orbits ; this inclination is therefore 
very nearly constant, as P'P^' is small and constant and very 
near Z. 

85. To integrate the eqvMion for the longitude of the 
epoch. 

We have (Art. 61) 

(JLt 

-f A^ee' cos (tsr — -bj') + A^ tan i tan i* cos (fl — H'). 
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Now from the fonmilse of Art. 72, we obtain 

e' = Jlf,' + M* + 2M^M^ cos {(^r, _ gr,U + 7^ - 7,}, 
e''=M^'+M,''+2M^M:cos {(g, -ff^t + y,- 7.}. 
ge' cos («r - w') = M^M^' + MJd; 

+ {M^M; +M^M^') cos {{g- -g,)t+y,- 7J. 

In like manner, &om the formulae of Art. 79, 

tan' i = N* + JV/ + 2iV;iV; cOs {hj; + S,- S j, 
tan" i'=N'^''+ N* + 2if/iV; cos [hj + S, - 8 J, 
tan i tan i' cos (n - il') = -ZV^, AY+ -ZV,* 

+ a; (i^, +i^/) cos {^ + 8, - n 

If these values be substituted in the expression for -5- , 
it takes the form 
rfe 

where Bn, B^^ and B^ denote certain constants. Integrating, 
we have 

6 = 6, + 5»< + — -L- sin {(5f, - flr,) « + 7, - 7,} 

+ |i sin {A.e - 8. - 8,}. 

We may omit the term Bnt, if we consider it as furnishing 
a correction on the mean motion n, which thus becomes 
(1 + B) n. With this understanding 

B B 

€ = €^+ -—^ sin \{ff^ ^g^t^ 7,.- 7J + T-' sin {h^t + 8, -SJ. 

y 1 "" if 9 '''1 

If this expression be developed, we may again omit the 
term involving the ferst power of ^, and consider it as affording 



90 PLANfiTAET THEORY. 

a farther correction to the mean motion*. Thus we shall 
obtain a series of the form 

86. In the Theory of the Planets this inequality is insen- 
sible, but in that of the Moon it amounts to upwards of 10 
seconds in a century, forming what is termed the secular ac- 
celeration of the Moon's mean motion. Thus it appears that 
this inequality does not, as its name would seem to imply, 
contradict the general theorem of the invariability of the mean 
motions, since it is due to a variation, not of the mean motion, 
(as we have employed the term in the preceding pages,) but 
of the epoch. If, however, as in the Lunar Theory, the epoch 
be omitted, any variation in the mean longitude will of neces- 
sity be thrown upon the mean motion ; only in this case, n 
will not be given by the equation n^a* =iu 

87. We have hitherto supposed the planetary motions to 
be referred to a fixed plane, but have left the particular plane 
undetermined. In practice it is usual to take the position of 
the ecliptic at some given epoch, as for instance the year 
1800 ; but since it is to the true ecliptic that astronomers 
refer the celestial motions, we will now obtain formulae for 
determiniDg relatively to the plane of the Earth's orbit, the 
position of that of any other planet. 

Let then m, m' denote the masses of the Earth and the 
planet considered, and suppose the orbits of m and m' but 
little inclined to each other and to the fixed plane of reference. 
Let X, X denote the latitudes of points in these orbits corre- 
sponding to the same longitude 6^\ then (see fig. to Art. 13) 

tan X = tan i sin (0^ — fl), tan \' = tan % sin {0^ — Of). 

* The advantage of thus diBposing of these terms arises from the fact 
that the mean motion, aci determined by observation, is the complete coeffi- 
cient of t in the expression of the mean longitnde. 
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Now since i and V are very small, we may replace tan X, 
tan X by X, X' respectively : tbvm 

X' — X = tan I sin (0^ — fl') — tan z sin (0^ — fl) 

= (tan i' cos fl' — tan i cos fl) sin ^^ 

— (tan i' sin fl' — tan i sin fl) cos 0^ , 
or, with the notation of Art. 79, 

V - X = ii'-q) sin ^i- (p' -p) cos 0^ (1). 

Now let 7 denote the inclination, v the longitude of the 
node of the orbit of mf relatively to that of m ; then approxi- 
mately 

X' — X = tan 7 sin (0^ — v) 

= tan 7 cos 1/ sin ^j — tan 7 sin 1/ cos ^j (2). 

Hence equating coefficients of sin 0^ and cos^^ in equations 
(1) and (2), 

tan 7 cos 1/ = g' — q, tan 7 sin v =p'— p ; 
whence tan' 7 = (p' — p)* + (? — q)\ 

and tan v = ^-? — ^ . 

q -2 

These expressions determine the position of the orbit of 
m! relatively to that of m, when the values of p, p', q, and q 
are known. Differentiating them, and neglecting small quan- 
tities of orders higher than the second, we obtain 

dy _ (dp' dp\ . (dq dq\ 

dv ^/d]^ dp\ cos V Idq dq\ sin v 
dt^\dt dtj tan7 \dt dt) tan 7 ' 

If the values of -^ , -^, &c. be substituted, these equa- 
tions give the variations of 7 and v. 
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88. It has been proposed to employ the invariable plane 
of the Solar System as a plane of reference ; and since it 
remains absolutely fixed, it would afford a means of com- 
paring together observations separated by long intervals. 
But the want of a convenient fixed line upon it, from which 
to measure longitudes*, stands in the way of its practical 
application ; so that it must be looked upon, at all events for 
the present, rather as theoretically interesting than as avail- 
able for astronomical purposes. 

* See Pont6coulant's Systlme du Monde, tome i. p. 469 (2ncl edit.). 

Note, Gauss in a memoir (Determinatio Attractionis^ dtc.) contributed to 
the Transactions of the Royal Society of Science of Gottingen, 1818, has 
stated without proof the following theorem. The secular variations which 
the elements of a planetary orbit experience from the perturbation produced 
by another planet are independent of the position of this planet in its orbit, 
and would be the same whether the disturbing planet nioved in its elliptic 
orbit according to £epler*s Laws, or whether its mass be conceived to be 
equably distributed throughout the orbit in such manner that equal parts of 
the mass are now assigned to such parts of the orbit as were described in 
equal intervals of time, provided that the times of revolution of the disturbed 
and disturbing planets are not commensurable. Gauss has also shewn how 
to reduce the expressions for the attraction of such a mass-orbit at any 
assigned mass-point not on the orbit to calculable forms of elliptic inte- 
grals, [a. f.] 



CHAPTER VL 



PERIODIC VARIATIONS OF THE ELEMENTS OF THE ORBIT. 

89. We come now to consider the variations produced 
by the periodical terms of R, These are called Periodical 
Variations, as opposed to the Secular Variations produced by 
the non-periodic terms. We have seen, indeed, that the 
latter are for the most part periodical in form, but in the 
Planetary Theory, the term Periodical Variations is restricted 
to those we are about to consider in the present Chapter. 

90. We have seen (Art. 46) that the general type of a 
periodical term is Pcos {{pn ± qri) t +Q}, where P is a func- 
tion of a, c, i ; and Q is a function of tsr, €, fl. Now such a 

.,, 1 • •! . • dR dR T dR 1 

term will produce a similar term "^ ;7~ , ~t- , and -jv-; but 

a term of the form P^m{{pn ±qri)t-\'Q} in -^, -7-, 

and -^, If, then, these be substituted in the equations 
of Art. 39, they will take the forms 

^ = P^An{{jpn±qn')t+Q}, 
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^ = P, COS {{pn ± qn') t + Qi 

de 

^ = P^ cos {(pn ± qn) t + Q}, 

fin 

^ = Pj cos {{pn ± qn') t + Q}, 

di 

^= P^sin {(pn ± qn") t-¥ Q}, 

d*t: 

^ = Pj8m{(pn±qn')t + Ql 

where P^, P^, &c. are functions of the elements of the dis- 
turbed and disturbing planets, and involve the first power of 
the disturbing mass. 

91. In integrating these equations, we may in general 
consider the elements which enter in the right-hand members 
as constant and equal to their values at the epoch from 
which the time is reckoned*. 

Let then a, e, «r, &c. denote the values of the elements 
at epoch, Sa, Se, Sw, &c. their periodical variations after an 
interval t : then integrating the above equations, and omit- 
ting the constant terms, we have 



pn± qn 
pn±qn 



Sa = i^fx>a{(pn±qnr)t + Q], 



Se = ^;r^' cos {(pn ± qn') t + Q}, 



p 

^ pn ±qn' ^^ ^^^ * *"'^ * "'' ^^' 

* This is eqniyalexit to neglecting ihe square of the disturbing force : see 
Art. 95. 
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Si ^/C08{(p»±2n')< + Q}, 



^^^ ~ (pn ±qnr ^'"^ ^^^ - ^^ ^ "^ ^^- 

Hence it appears that the variations produced by the 
periodical terms of JB are all periodical in form. 

92. It will be seen that all the expressions of the last 
Article involve the divisor pn ± qn\ while Sf involves the 
divisor (pn ± quf. If then it should happen that either 
pn + qn or pi - qii' is very small, a term in JB containing 
(pn ± qn') t in its argument, though of a high order, may 
have a sensible eflfect on the elements of the orbit. Now 
since p and q are either positive integers or zero, pn + qn' 
cannot be small unless n and n' are small, a case which does 
not occur with any of the planets : but we have instances 
in which pn - qn is small*. 

Since the period of such inequalities is very great (being 

,] , they are called inequalities of long period, or long 

inequalities, 

93. To select such terms in R as will produce the prin- 
cipal inequalities of long period. 

* If in any case the mean motions of two planets were exactly commen- 
surable and in the ratio of i? to q, the corresponding term of B, as we have 
already remarked (Art. 48), would cease to be periodical and would form 
a part of F^ but no instance of this occurs among the planets. 
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We have seen that the dimension of the principal part 
of the coefficient of a term containing (pn ^^ qnT) t in its argu- 
ment is jp *- g' (Art. 50) ; hence if we can find two integers 
p and q nearly in t1ie ratio of n to n, and having a small 
difference, the corresponding term of R will produce an im- 
portant long inequality in the elements of each planet. 

In the case of Jupiter and Saturn w : w' :: 5 : 2 nearly, 
and 5—2 = 3; hence there is a long inequality arising from 
a term in R of the form Pcos{(2n — 5w') *+ Q}, the prin- 
cipal part of P being of the third order. This inequality is 
interesting in an historical point of view, having long baffled 
the labours of mathematicians and appeared inexplicable on 
the hypothesis of gravitation. It was at last successfully 
explained by Laplace, 

For the Earth and Venus, n : n' iiS^ilS nearly, so that 
there is a long inequality arising from a term in R of the 
fifth order. The discovery of this inequality is due to the 
Astronomer RoyaL 

Finally, in the case of Neptune and Uranus, n :nf ::1 : 2 
nearly, hence there is a long inequality arising from a term 
in R which is of the first order. 

94. Between corresponding terms of the long inequali- 
ties in the mean motions of two planets, arising from the 
near commensurability of n and n, there is a simple approx- 
imate relation. 

Let m, m' be the masses of the two planets, R, R' their 
disturbing functions : then by Art. 8, considering only the 
mutual action of m and m\ we have 



?w' m' 



r 
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^ We shaU distinguish the first and second terms of R and 
R' as the symmetrical and unsymm^rical parts respectively, 
since the co-ordinates of m and m' are involved symmetri- 
cally in the former but not in the latter. 

Since, then, the symmetrical parts of R and R' differ only 
in having m and m' interchanged, if 

m'if cos {(pw *- qn) t-\'Q] 

be any term in the symmetrical part of iJ, that of iJ' will con- 
tain the term 

^if cos {(pn — qn') t + Q}, 

Confining our attention to these terms, we have (Art. 39) 

rf"f 3n'a dR Snad, (R) 
cf^' jjL de fh dt 

= ^^^ m'M^\n{{pn — qn)t+Q]y 
therefore Sf = - -^-^ 7 tts sin [{pn — an') < + Ql. 

Similarly, 8^= ^ (^^^ «^ i(^^ " «^') ^ + «}• 

Sf m'filn^ap m'fina 

Hence srr? = 72~^ ~ r~f 

oj m^in^aq m/juna 

approximately, since qn' is nearly equal to pn ; therefore 

8g ^ m' V(mV) 
Sf mVO^a) ' 

or, since /a' differs from /a by a quantity of the order of the 
c. P. T. 7 
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disturbing force, the square of which we are neglecting, 
we have to the first order^ 

the required relation. 

The same relation is also approximately true in the case 
of terms arising from the unsymmetrical parts of jB and jB'*. 
For denoting these by iJ^ and ii/ respectively, we have 

„ , ( x' y' ^\ 



i?/=-m(4 + /J + /^). 



Now the equations of motion of the planet m' referred to 
rectangular axes are 

dt' "^ t"" " dx' '^^*' 

and since, by the principles of the method of the variation of 
elements, the /orw of solution of these equations is the same 
as it would be if JS' were zero, it follows that, if the differ- 
ential coefficients be taken a« if the elements were constant, 
we shall have 

r» / df ' r"" iif df ' r'*"" /t' df ' 
and therefore, with this understanding, that 

^' (I'V df ^y de ^' dfj- 

* For the demonstration of this we are mainly indebted to The Theory of 
the Long Inequality of Uranm and Neptune: an essay which obtained the 
Adams Prize for the year 1850. By E. Pierson, M.A. 
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Similarly. R^=-[a>-^^y - + , -^ 

Now any term in R^ containing (pn — qn) t in its argu- 
ment can arise only from the combination of terms in x, y, 

and z, containing pnt in their arguments with terms in -^-^- , 

dW d^z 

-^> and -175- containing qrit Suppose, then, x and a?' when 

developed in terms of t and the elements to contain respec- 
tively the terms 

L cos {^pnt + I), U cos {qrit + V). 

d^od 
Hence the product x -r^ will contain the term 

- \ LL'ifn'^ cos [{pn - g/i') « + Z - r}, 

d^x 
and the product oi -^ the term 

— ^ LUpW cos {{pn — qn) t + l—V]'. 

the coeflBcients are in the ratio ^n^ to jpV. Similarly, the 

d^V d^z 

coefficients of the same cosine in y ——■ and z -^ are to 

d^V d^z 

those in y' -7^ and z' -^ in the same ratio. 



Hence if -7- Mq^n^ cos {(pn — g^n") t-^Q} 

be any term in JJ^, then JJ/ will contain the term 

m 



ilfpV cos {{pn — gyi') ^ + Q], 



7—2 



L 
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Confining our attention to these terms, we have 

T- St m'aq m'na , 

Hence k^, = , = j-r nearly, 

Sf majp mna ^ 

m Aja ' 

the square of the disturbing force being neglected. 

By means of this relation, when one of the long inequali- 
ties is known, the other may be calculated : it may also be 
used as a formula of verification. 

95. We have remarked that in integrating the equations 
of Art. 90, we may in general consider the elements which 
enter in the right-hand members as constant and equal to 
their values at the epoch from which the time is measured. 
In the case, however, of inequalities whose periods are very 
long, the secular variations of the elements in the interval 
produce a sensible effect. In order to take account of these, 
we may integrate our equations by parts, considering the 
elements variable ; and then substitute their values as calcu- 
lated by the method of Art. 62. For example, consider the 
equation 

^^ => P sin {{pn - qn") t + Q} 

= P sin X, suppose. 

Integrating by parts, and remembering that n is constant 
with regard to secular variations, we have* 
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-T2 = 7 COS X + ; T-j -J- Sin \ 

^ 1 _ cPP 

/\8 



{pn-qny df 



cos X — ... ; 



therefore 8? = — 7 ^^5 sin X — -. fr-^ -r: cos X 

\'pti, — qp!) {pn-qny dt 

1 cTP . _ 

cos X + 7 TTi — rs-sm X + . •. 



(pn-q;ny dt ^^'^^ (pn-^qny df 



1 d^P . ^^ 

sinX + ... 



{pn—qviy^ di 
d^P 



— .,,> sinX 



^f P 3 d^P 

1 {pn — qny {pn — qn'f di^ 

{ 2 dP 1 

' dP d^P 

In this equation P, -^ , -^ , &c. are functions of the 

elements ; their values may be calculated by the formula3 of 

dP 
Art. 62. It may be noticed that P is of the first order, -j- of 

(?P 
the second, and -^ of the third of the disturbing force : 

dP 

for -J- , being found from P by differentiation, will involve 

the differential coefficients of the elements, which are them- 

cPP 
selves of the first order ; and similarly for -7-5- . 

96. Having now completed our account of the methods 
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of treating the secular and periodic variations of the elements 
of the orbit, we will say a few words on the distinction be- 
tween them. In the first place we may observe that the 
periodic variations involve the mean longitude of the dis- 
turbed and disturbing planets, and therefore depend chiefly 
upon the configuration of the planetary system. On the con- 
trary the secular variations depend solely upon the values of 
the elements. The latter class of variations take place with 
extreme slowness, so that if these only existed, a considerable 
time must elapse before the deviation of the planet from 
elliptic motion became appreciable. On the other hand, the 
periodic variations (such at least as are rapidly periodic) 
"are in their nature transient and temporary: they disappear 
in short periods, and leave no trace. The planet is tempo- 
rarily drawn from its orbit (its slowly varying orbit), but 
forthwith returns to it, to deviate presently as much the other 
way, while the varied orbit accommodates and adjusts itself 
to the average of these excursions on either side of it ; and 
thus continues to present, for a. succession of indefinite ages, 
a kind of medium picture of all that the planet has been 
doing in their lapse, in which the expression and character 
is preserved; but the individual features are merged and 
lost*." On this account it is convenient to suppose the 
planet to move in an ellipse, the elements of which are cor- 
rected for secular variations only, and to take account of the 
periodic variations by applying small corrections to the radius 
vector and longitude as calculated from the elliptic formulae. 

97. We will accoirdingly shew how by means of the 
periodic variations of the elements, the corresponding in- 
equalities in the radius vector and longitude may be calcu- 
lated. If we' take for our plane of reference the position of 

• Herschers Outlines of Astronomy, 10th edit. Art. 656. 
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the plane of the orbit of the disturbed planet at the epoch 
from which the time is reckoned, the inclination will be of 
the order of the disturbing iorce, and therefore, if we neglect 
the square of the latter, we may also neglect the square of 
the former. 

98. To calculate the periodic variations in radius vector. 

Let Sa, he, S«r, .&c. denote the periodic variations in 
a, c, -Bj, &a, and let hr be the corresponding variation in r ; 
then 

^ dr ^ dr ^ , dr ^ dr ^^ dr ^ 

in which the square of the disturbing force is neglected, 

since this would be introduced by the squares and products 

of ha, he, &c. The values of Sa, Se, &c/have been found: in 

dr dr 
Art. 91, those of -j- , -r- , &c. may be obtained from the 

equation (Art. 40) 

f 1 1 1 

r = a jl + 2 e*- ecos (f + € --or) - ^ ^ cos 2 (f + e - ut) - ... •. 

99. To calculate the periodic variations in longitude. 

These might be found in the same manner as the varia- 
tions in radius vector, but they may also be deduced from 
them: we proceed to obtain a formula for this purpose. 

We have ^ = ^' (Art 22), 

and ^-^o = f2-Xl^; 

therefore "5/ ~ "^ "^ ^"^ "" ^^^ *^ "^ ^^^^ ^^' ^^^ 

Jh 
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since (1 — cos i) -i- , being of the order of the cube of the 

disturbing force, may be neglected. 

Let Sr, BO, and hh be corresponding variations m r, 
and h ; then 

dt {r + ^)*'' 



de die 



or -^ + 



-^(^4*)('4J 



_h ih ihir 

neglecting the square of the disturbing force ; therefore 

di6_hh 2A8r 
dt f* 7^ ' 

which gives the variations in longitude. The value of Zh 
may be found from the formula 

dh_dR dR 
dt de d'os * 

For the periodic variations in latitude, we refer to Pont^- 
coulant's Syst^me du Monde, Tome i. p. 492. 

100. As an example of the processes of this Chapter, we 
will calculate the variations in radius vector and longitude 
due to the term mMe cos {(n — 2n) ^ + € — 26' + «r} in jB. 

Considering this term only, we have 

jB = m'Me cos {(n - 2nO « + e - 2€' + wj 
= m!Me cos \ suppose. 

Hence 



dR ,dM 

-r- = w -^— e COS A., 
da da 


dR , , -. ^ 

de 


dR ,,^ . ^ 


dR ,,^ . _ 

-— sr — mMe sm X, 


di^' 


do. "• 
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Substituting these in the fortnula& of Art. 59, and neg- 
lecting small quantities of orders higher than the first, we 
have 

da 2nc? ,,^ . ' 
at II 

de Tia i-mr • ^ 
— -=s — mirf sinX, 
at !>» 

dw na ,,-. _ 

e-7r= — mMcosX, 
at fi 

de 2na^ ,dM ^ 1 na ,,, 

— .^ _^ e cos X + H — emm cos \, 

at fi aa z fi 

zr^ = • mMe smX. 

air fi 

By integration we have 

^ 2mM na*e 

oa = ZT"/ cos A, 

/i n — 2w- 

-, tnM na 

oe = — TT"? cos A, 

fi n — zn 

^ 7n!M na . ^ 

coo* = 2n sin a, 

/i n — 2m 

' ^ /Im'M 2madM\ nae . ^ 

\2 /A /A dajn — zn 

^^^ SmfM n^ae . _ 

5f — ^ ^-^^sinA. 

fi (n — zn) 

Also 
r = ajl +se'^-^cos(5'+e — «•)— ^e*cos2(5'+€ — «^) — ...[; 



106 ^LANETAET THEORY. 

therefore, small quantities of orders higher than the first 
being neglected, 

^=l-eC0S(f+€-iar), 

civ 

^ = a {e — COS (5"+ € - ot) - e cos 2 (5"+ e — -bt)}, 

dr 

^ = — a {« sin (5" + e - ut) + e* sin 2 (f + € — isr)}, 

^=a6Sin(f+e--B5-), 

^= ae sm (f + € - tsr). 

XT r. dr ^ dr ^ dr ^ dr ^^ dr ^ 
da de acr d^ de 

2mM nd?e 

-, COS X 



/A w — 2n 

8 

7 sin X {sin (f + e — tsr) + e sin 2 (g" + e — w)}. 



Tn'Jf wa 



m'ilf na 



Since we are neglecting the square of the disturbing force, 
the elements in this equation may be considered as constant, 
and therefore nt written for ^: we have then, restoring to \ 
its value 

Sr = ^r-j cos 2 [ (n — w J < + e - € } 

fi n — zn ^^ ' 

+ ^r-, cos {(n-2wO^ + e-2€ +tsr} 

fi n — zn ^^ ' , 



H.w~ 
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H ^— , cos {(3n — 2n)^ + 3e-2e — tsrj, 

which is the variation in radius vector. 



101. To calculate the variation in longitude, we shall 
employ the equation 

dhe hh 2hZr 
dt " r" r* • 

-^T dh dR . dR e% f-kJT • % 
Now j- = ^- + J— = — 2mJfesinX, 
dt de dnr 

therefore Sh = — -^ri cos \ ; 

n—zn . 

^, „ SA 2w'Jlfc 

therefore -«- = -r^ t-tt cos X 

r' a'* (n — 2n ) 

-,cos{(n--2n)^+e — 2e +«rj. 



/Lt n— 2n 
Also — j- = — 5-{l + 3ecos(7i< + e--'BT)H- ...} 

= ST COS 2 {(« — W ) ^ + € - € } 

/A n — 2/1 ^^ ' ^ 

"^2"/r;r:^'^^'f(^"2^^*+"-^'+^J 

+ 2 7- iT^l^' *'**' t(3« - 2n ) « + 3€ - 2e - «rl. 

Hence by substitution 

die 2m M tfa ^.. ,. . , ^ 

^=— 7r-^r::2;^'^"'2{(n-n)e+6-e} 
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Sm'M iffae 



fi w — 2n 

5m'M n^ae 
fi w — 2n' 

By integrati6n 



7COS {(n-2n')< + e-2€' + Br} 



7 cos {(3n - 2»') « + 3€ - 26' - w}. 



g^= ^ TV sin 2 {(n-n') <+€-€'} 

ft (n — 2n)(n — n) ^^ "^ ' 

7 Q-Txa Sin {(n - 2n') « + € - 2e' + ^} 

- A* (»-2n-)(3»-2nf ^°t(3>>-2«)t + 3e-2e-H 
which is the variation in longitude. 

In the case of Uranus and Neptune, since n : ti nearly 
as 2 : 1, the term we have been considering is important in 
the theory of the long inequality. 



CHAPTER riL 

DIRECT METHOD OF CALCULATING THE INEQUALITIES IN 
RADIUS VECTOR^ LONGITUDE, AND LATITUDE. 

102. In the calculation of the planetary inequalities, we 
have hitherto employed exclusively the method of the Varia- 
tion of Elements, but there is another method of solving the 
problem, which demands our attention. It consists in obtain- 
ing equations for calculating the inequalities in radius vector, 
longitude, and latitude directly from the equations of motion. 
This method is indeed the simplest to employ in the case of 
periodic variations of short period, that of the preceding 
Chapters being the most convenient for the calculation of 
secular variations and long inequalities. For, since these 
latter take place with extreme slowness, the elliptic ele- 
ments, when once corrected for them, continue for a consi- 
derable period to represent the actual motion ; while, in the 
former case, the values of the elements change rapidly, and 
the motion cannot for long be represented by the same 
ellipse. We proceed, then, to the direct method of calcu- 
lation*. 

103. If rj, 6^, and z denote the projected radius vector, 
longitude, and distance from the plane of reference, of the 
planet, we have (see Art. 9) the equations of motion 

* The two methods are sometunes distinguished as those of Lagrange 
and Laplace, but in the M€canique Cileste we find both employed. 



110 PLANETARY THEORY. 

d\ fdd\* fir, , dS 



df '\dtj r» ' dr^ 

( »£^A_^ 
V' dtj~de. 



d'z _ /tt.8 dR 
de~~l''^dl ' 

If we take for the fixed plane of reference the position of 
the plane of the orbit at the epoch from which the time is 
reckoned, the inclination (as we have remarked in Art. 97) 
will be the order of the disturbing force, the square of which 
will be neglected. Now it will be seen on referring to Art. 
42, that r^ and 0^ differ from r and by quantities depending 
upon the square of the inclination : hence in the above equa- 
tions, we may replace r, and 0^ by r and respectively. Also 
if \ denote the latitude of the planet, we have 

« = r sin \. 

Hence our equations of motion become 

<fr fd0\^ fi dR 



d 
dt 



(^,d0\dR 
d'(rsinX) _ /^ siTi\4-— n^ 



df r" dz 

104. As a first approximation, let values of r, and X 
be obtained from these equations by neglecting the disturbing 
force, and let r H- Sr, 5 + S5, \ + S\ denote the true values of 
these co-ordinates; then Sr, h0 and SX will be very small 
quantities, of the order of the disturbing force: they are 
termed the perturbations in radius vector, longitude, and 
latitude. We proceed to investigate equations by means of 
which these quantities may be determined. 
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105. To obtain the equation for the perturbation in radius 
vector. 

From equations (1) and (2) of Art. 108, we obtain 
UJ U/ r^J\drdt^dddt)^ 



-^^/ 



^dt+C (4). 



Multiply (1) by r and add to it (4) : thus 



r 

I 

I 



j-(i)"-?-f-^/^*-''. 



If the disturbing force be neglected, this equation be- 
comes 

^-i*-^" w- 

Let a value of r be obtained from this equation, and let 
r + Sr denote the true radius vector : then if we agree to 
neglect the square of the disturbing force in our next ap- 
proximation, it will be suflScient to write r-{-hr for r in those 
terms of (5) which do not involve the disturbing force : also 
since hr is itself of the order of the disturbing force, its square 
may be neglected. Hence from (5) 

therefore - ,;» + 2 — \^ « -J- — ^ Sr-h 2r -p 
dt dv r r\ . dr 
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hence by (6), 

which is the equation for the perturbation in radius vector. 
We may express the right-hand member in a more convenient 
form, for since (Art. 42) 

r = a(l + w), 

dR __ dB ^^ _ /■■ , \ ^R . 
da " dr da"^ ^ dr ^ 

, ^ dR dR 

therefore . ^ j— = o^ 3- ; 

dr da 

d(R) dR 
Hence our equation becomes 

106. To obtain the eqvutionfor the perturbation in longi- 
tude. 

We have from equation (1) of Art. 103, 



/^'^^IdV u _1 



dR 



\dt J r d^ r* r dr ^ 

As before, let a value of 6 be obtained from this equation, 
the disturbing force being neglected, and let d + hd denote 
the true value of the longitude: then writing r-\-Zr for r 
and 6 -¥^6 for 5, and neglecting the square of the disturbing 
force, we have 

dddW^l^ ^(^_Sfi^ _ldR 
dt dt "r de^ r^ de r* r dr' 

_ . ^dd 7 dR dR, 
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therefore 2&-jr-=r->ia — Sr ^-i — ^or-^a-j- 

dt dv dv r da 

d / dSr ^ dr\ 3/jl ^ dR 

This equation will become integrable if we eliminate the 
term — -^ . rBr by means of the equation for the perturbation 
in radius vector. We have from that equa4}ion 

th^efore by addition. 



n [dR , 

therefore m0 = Z^ + r^-^^-^{^ dt 

dt dt dt J da 



the arbitrary constant being considered as included in the 
sign of integration ; therefore 

^-*'-^-^4'-*'/S*-''"//f'*'> 



or 



^'''-^-^%-Ht^-^ii^^- 



which determines the perturbation in longitude, when that in 
radius vector is known. 

c. P. T. 8 
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107. To obtain the equation for the perturbation in 
latitude. 

From equation (3) of Art. 103, we have 

(^(rsmX) /A (r sin X) '_ c£R 
df ^ ? ■&• 

Since the plane of reference is supposed to coincide with 
the position of the plane of the orbit at the epoch from which 
the time is reckoned, we have at the epoch \ == : hence, 
denoting by SX the latitude at time t, our equation becomes 

which is similar in form to the equation for the perturbation 
in radius vector, 

108. To integrate the equation for the perturbation in 
radius vector. 

The equation is (Art. 105) 

d^irSr) , A*/ 5S . dR^^ {dR .. 

Let us consider a term in 22 of the form 

P cos {{pn — qn*) t + Q}, 

where P is a function of the mean distances, excentricities, 
and inclinations, and Q of the longitudes of the perihelia, 
nodes, and epochs : then uniting this term with the non- 
periodic part of JB, which we have denoted by F^ we have 

^ = P + Pcos {(2w - gnO ^ + Q} ; 
therefore ^f = ^+^cos{(pr.-3ri')« + «l. 

g=-.pgsin{(pn-yn')< + (2}. 
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since F does not contain e ; therefore 



"/ 



dR ,. de 



where ^ is an arbitrary constant. It may appear superflu- 
ous to introduce this quantity, since an arbitrary constant C 
has ahready been added in Art. 105 ; but its introduction is 
merely equivalent to a small change in the value of C, and 
any alteration which does not interfere with the first ap- 
proximation, obtained by neglecting the disturbing force, is 
of course permissible. Thus ^r is a purely arbitrary quantity 
which might have been omitted, but which we shall find it 
convenient to retain, leaving its value to be assigned here- 
after. 

„ dR ^ ^ [dR ., ^ , ^ dF 

Hence a ^ + 2n I ^ a^ = zmg + « j~ 



+ 



.p ^-^f I 



a-7- + — r (cos [{pn'-qn)t + Q} 
= 2m g + a^ + jPj cos {(pn - jn') t -t- Q], 



Mppose, wkere -P. = « "3^ + JJTJJ • 
Again (see Art. 13), 
r = a]l + 2^-"^cos(n^ + 6-«T) 



1 
— ^e'cos2(n^ + 6 — -cr) — ... 



therefore 4 == 4 ll + 3e cos (w^ + €-tsr) + ...k 

r a [ ) 

= w* {1 + Se cos (ni + 6 — «r) + ...}, 
since nV « fi. 



8—2 
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Hence, by substitution, the equation for the perturbation 
in radius vector becomes 

^^^ + n".rSr = 2m> + a^+ P, cos (Ojn-gwO^+Ql 

— n' rSr {Se cos (n^ + e — cr) + ,..}. 

109. This equation must be solved by successive ap- 
proximation, as in the Lunar Theory. By omitting all 
small quantities, we obtain a first approximation to the value 
of rBr ; this being substituted in the second member, and 
small quantities of orders higher than the first neglected, we 
obtain a second approximation, which will be correct to the 
first order. In like manner, a third, and higher approxima- 
tions may be obtained. 

On referring to Art. 59, it will be seen that small quanti- 
ties of the second order being neglected, 

F= I m'C, ; 

hence, neglecting all small quantities, the equation of the 
preceding Article becomes 

The integral of this equation is 

+ • 2 — 7 — * K2 COS {{pn - 071") t + Q} 

+ ^cos(7i^— B), 
where A and B are arbitrary constants. Since, if all small 
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quantities be neglected, r = a, we have as a first approxima- 
tion, 

A 
H — cos {at — B). 

110. We may, however, omit the last term: for, con- 
sidering this only, the radius vector of the planet becomes 

a -jl — e cos (n< + € — -Bj) + -jcos (n^ — £)+ ... j- 
= a [1 — {e cos (e — w) — 5 cos B\ cos tit 

Cb 

A 
+ {e sin (e — w) + -^ sin B] sin w^ + ...] 

= a {1 — gj cos {nt + e — -nrj 4- ...}, 

-4 

if ^j cos (e — isr J = e cos (e — -Bj) ^ cos 5, 

A 

e^ sin (e — wj = e sin (e — w) 4- — j sin £, 

from which e^ and w^ may be determined. 

Now since the ellipse upon which our approximations are 
based, has been obtained by neglecting the disturbing force, 
we may in the elliptic formulae replace e and -cr by e^ and isr^ 
respectively, since they differ by quantities of the order of the 
disturbing force. If this be done, pur first approximation 
becomes 

p 

j-^. cos {[pn — qn') t + Q}. 



a [ri^ — [pn — qnf] 
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111. la order to obtain a secoad approximation, this 
value must be substituted for Sr in the right-hand member of 
the equation of Art. 108. Also since the square of the dis- 
turbing force is neglected, we may write Cj and w^ for e and w 
in this equation. We will write for brevity 

Br = L + P^ cos {(pn — jn') t + Q}. 

Substituting this in the equation of Art. 108, and omit- 
ting those terms which have produced the first approxima- 
tion*, we have 

,w +7r.rSr 

= — 3w"a«j cos {nt + € — •cTj) [L + P^ cos {{pn — qn') t + Q}] 
= — ^ri^ae^ L cos,(n^ + e — •orj ' 

3 

- ^n'a^^P, cos [{{p + 1) n - qn']t'{- ^ + €- tsrj 

S 

— 2 w'o^i P J cos [{ (p - 1 ) w - jn'} ^ + Q - € + «r J . 

112. On the form of this equation, we have an import- 
ant remark to make. In consequence of the term 

— Sn'o^j L cos (w^ + € — -BTj), 

its integral will contain the term 

3 

^ ae^ L cos (n^ -h e — Wj), 

Thus we are met by a difficulty: our equations have been 
formed on the hypothesis that the square of Sr is small 
enough to be omitted, whereas here, we have a term capable 
of indefinite increase. This term, then, if retained, would 
ultimately vitiate the whole approximation. The difficulty 
might, as in the Lunar Theory, be obviated by writing 

* These tenns are omitted for the sake of brevity : in order, therefore, to 
obtain the complete second approximation, we mast add to the integral of the 
above equation the result of the first approximation. 
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cn for n in the elliptic formul^^ which amounts to supposing 
the perihelion to be in motion. Its motion is however better 
found by the method of the variation of elements. Indeed it 
may be shewn that such terms as those we are considering 
lead to the formulae which have already be^n obtained, for the 
secular variation of the elements*. We may accordingly 
altogether neglect such terms if we suppose the elements of 
the ellipse on which our approximations are based to have 
been previously corrected for their secular variations. 

113. With this understanding, the complete integral of 
the equation of Art. Ill will be 

5. 3 n^ae, P^ 
fCr = — -^ ' 



cos[{{p + l) n- qn'} t + Q + e-'srj 

3 n^ae, P, 

2n!''-{(p'-l)n-qnf] 

cos[{(|)-l)n-gn']« + Q-€ + t!rJ 
+ A cos {nt -^ B). 

If this be added to the result of the first approximation, 
we obtain for a second approximation 



(pn — qny 

cos {(j>n — q9i) t + Q] 



n^ae^ P, 



2n^^{(p + l)n^qnf 

cos [{(p + 1) yi - qnlt + Q + € - t!rj 

3 n^ae, P, 

2n^^{{p-l)n-qnY 

cos [{(jp — 1) n — qn'] t+Q — e + xsr^] 

+ -4 cos {nt — B\ 

* It is thus tliat the Secular Variationfi are first obtain ed in the Micanique 
Celeste, See Pont^oulant's Systeme da Mande, Supplement au Livre IL 



^mitmutmmifmmmmmm^m^mmi^mmmmtmmmm^m^^^^m^c^m~^^^9mm^^^mt^^^^^^^^^\r ■ ■ ■■ 
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The arbitrary constante may be disposed of as in the first 
approximation. In order to obtedn a second approximation 
to the value of 8r, it is only necessary to multiply the right- 
hand member of the above equation by 



-• {I + ^1 cos (nt+ € — ts-,)}, 

Cb 



neglecting e^\ 



114. To calculate the perturbations in longitude. 
We have (Art. 106) 

Taking for simplicity, the first approximation to the value 
of rSr, which has been obtained by neglecting the first power 
of the excentricity, we have 

eoB{{pn--q'nr)t+ Q); 
therefore 

^ d . rSr 2P, (pn- an*) . ,, ,^ ^ r\^ 

also, neglecting the first power of the excentricity, 

I 
/ 

As in Art. 108, writing^ ^rn'O^ for Fy we have 
dR 1 ,dC. dP „ .^ -, 
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dR ,. ^ ffdR 



therefore — 2a I -r- d^ — 3n 1 1 



de 

, da 



dt 



.9 



pa — qn^ da jpn^qn f 



pn — qn 
where /ii^ an arbitrary constant. 

Hence by substitution^ we have 



8in{(pn-qn')t+ Q], 



hSe =/- (m'a ^ + 3wV) t 



2a f 3.pf 
da de 



[jpM — jw {pn—qnY 

115. This expression is open to the objection of contain- 
ing a term proportional to the time, which being capable of 
indefinite increase, would nltimately vitiate the whole ap- 
proximation. Here, then, we see the advantage of having 
a quantity g which may be determined at pleasure : we will 
so determine it that the objectionable term shall vanish. 
This condition gives 

1 dc. 
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We may also omit the constant/, and consider it as con- 
tained in the epoch. We have then, writing for h its value 
na^ V(l — e'), and neglecting e", 

A^af 2nPf 
^yj 1 < aa , de 

na \pn — qn {pn — qn) 



2R 
n 



?§|C^)^n|(.« -«. + «). 



116. Before proceeding to obtain the perturbations in 
latitude, yre will make a few remarks on the forms of the ex- 
pressions for Sr and B0. If we confine ourselves to the results 
of the first approximation, it will be seen on substituting the 
value of Pj, that pn — qn and n' — (pii — qny occur as divi- 
sors, and that the expression for Bd contains besides, the 
divisor (pn — qnTf. The second of these may be written 

{(1 -p) n + qn'\ {(1 + p)n - qn]. 

If then either 

(i) pn — qn', (ii) (1 — p) w 4- qn\ or (iii) {l+p)n — qn\ 

be very small, the corresponding terms in Sr and B0, though 
of a high order, may yet be sensible. This is especially the 
case with the first, since as we have remarked, its square 
occurs in the expression for B0. These are instances of what 
in the preceding Chapter have been characterised as long 
inequalities. 

The period of the term P cos {(pn — qn") t + Q], which has 
given rise to these inequalities, is 

2v 
pn --qn'' 

in the case of (i), this is very large, and in that of (ii) or (iii) 

27r 
it is very nearly equal to — , since pn — qn' is nearly equal 
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to + n. Hence it appears that terms in 2J wliose period is 
either very large, or nearly equal to that of the planet, may 
give rise to important inequalities in the radius vector and 
longitude. Their actual importance will of course depend in 
part upon the order of the principal part of P with respect to 
the excentricities and inclinations, ie. (see Art. 50) upon 

117. To integrate the equation for the perturbation in 
latitude. 

The equation is (Art. 107) 

the position of the plane of the orbit of the disturbed planet 
at the epoch being taken for the fixed plane of reference. 

Differentiating the expression for B in Art. 44, with respect 
to z, we obtain 

- m(z — z') (iD^ + Dj cos + ... + 2?^ coshj} + ...) 7^ , 

a 

or, putting z equal to 0, and substituting 

a' tan i' sin (nt + c' — fl') 
for z' (see Art. 42), 

iWtSini^sm(n't + €-a')i^D^^\ + D^cos<l>+..X. 

This expression, after reduction, consists of terms of the 
form 

Psin{(pi-57i')« + Q}, 

where p and q are positive integers, and either may be zero. 
Considering one such term, our equation becomes 



m 
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Now as in Art. 108, 

^ = 7^' {1 + 3e cos (n^ + e — -bt) + , . .} ; 

hence, neglecting the product eS\, we have for a first ap- 
proximation 

The integral of this equation is 

p 

rB\ = -o — 7 :^ sin {(pn — qn") ^ + §} +-4 cos (nt — B). 

If instead of taking for the fixed plane of reference, the 
plane of the orbit of the disturbed planet at epoch, we take a 
plane slightly inclined to this, we may omit the arbitrary 
term. For, denoting the planet's latitude with respect to this 
plane by X, we have approximately 

\ = tan z sin (n« + e — O), 

and it may be shewn as in Art. 110, that omitting the term 
in question is only equivalent to changing slightly the values 
of i and il. 



CHAPTER nil. 

ON THE EFFECTS WHICH A RESISTING MEDIUM WOULD 
PRODUCE IN THE MOTIONS OF THE PLANETS. 

118. In the preceding Chapters, we have supposed the 
planetary motions to take place in free space, and the results 
of calculations based upon this hypothesis manifest a very 
close agxeemeni with observation. There is, however, a 
remarkable circumstance connected with Encke's comet which 
seems to indicate the possibility of the existence of a very rare 
medium, too rare iniieed to cause any sensible resistance to 
the motions of the planets, but which, as we shall presently 
see, may yet influence the motions of comets, in consequence 
of ihe extreme smallness of the masses of these bodies. It 
has been observed that the comet above referred to (which 
describes an elliptic orbit in a period of about SJ years,) has 
since its appearance in 1786, been moving round the Sun 
with an increasing mean motion. Encke attributes this to 
the resistance of a medium pervading space. We shall 
therefore proceed to examine the effects which such a medium 
would produce upon the elements of a planet's orbit, assuming, 
in accordance with the usual theory, that the resistance 
varies as the product of the density of the medium and the 
square of the velocity of the planet. We shall neglect, in 
the present investigation, all forces except the Sun's attrac- 
tion and the resistance of the medium; consequently the 
planet may be supposed to move wholly in one plane. 
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119. Let r, ^ be the radius vector and longitude of the 
planet, 8 the length of an arc of its actual orbit measured from 
some fixed point to its position at time t, and p the density of 
the medium. Then if ft be a constant, we may represent the 

resistance on the planet by kp [-f,)> aiid the equations of 

motion will be 

de \dt) ~ i> "PKdt) da' 

1 d (,dd\_ . [dsV dd 
rdt\dt) ''P\dt)''d^- 

tJff 
If r" -IT = A, these may be writtea 

dC "^[dtj r* "f'dtdt ^'' 

l(-i)-^l- P)- 

These equations are the same in form with those of 
Art. 20, and may be treated in a similar manner, -^kp-r: -f. 

taking the place of -r-, and ^kph-r: that of -t^, (See 
Art. 24.) We have firom equation (2) 

dh 7 -L ds 

di=-^p^dt' 

120. To obtain a formula for calculating the mean dis- 
tance. 

We might proceed as in Art. 26, but we shall here employ 
the method of Art. 26. We have 



(P« fidr J /efoy 
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and by a known formula of elliptic motion 

\dt) '^ T a 
Differentiating the latter, we obtain 

ds 
and multiplying the former by 2 -^, 



therefore 



^dtde l^dt "^^PydtJ' 



or 



a^dt 
da 



2hpa^ fds\* 
fi [dtj' 



121, To obtain a formula for calculating the excentricity. 
We have, as in Art. 27, 



§-{%h^r-'^ (^)- 



dr ds 
Differentiating as if r were constant, and writing "~ ^-j1 j^ 

d»r 



^''dr 



fj?ede_ , fdr\*d8 (fh /A/l fi\ ft'e*) dh 
h'dt "'' \dt) dt "^ l\r h) [r ^hV **" A» J di 



de _ 2kp ds 
dt" e dt 
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Now from equation (3) 

[dtj "*"r*~ A" ■*" r A' 

2A?p (1 - ^) fa \ds 
e U" /dr 

This result may also be obtained by diflferentiating the 
formula A' = fui (1 — e'), and substituting the expressions for 

-T- and j7 , as in Art.. 28. 
d^ dt 

122. To obtain aformvlafor calcvlating ike longitude of 
perihelion. 

We have, as in Art. 29, 

|cot(^-.) = ^-£ (4). 

Differentiating as if r and were constant, and writing 
y drds J. ^r 

or since -=.-kph^^, 
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but from equation (4) of Art. 22, 

dr ,^ \ H^ 

^cosec(5-isr) = ^; 

therefore -^ = ^ sin (5 - «r) -=? . 

a^ 6 ^ ' dt 

123. To o6^ai» a formula for calcvlating the longitude 
of the epoch. 

We have (see Art. 13), 

Differentiating, the elements being considered variable, 

d& _dm \ fdf\ d (n< + € - "sr) d/*de^ 
dt "" di n. \c?^ / dt dedt^ 

and differentiating as if the elements were invariable, 

d0_/df\ 
dt^UtJ* 

Equating the two values of -i^ , we have 

dm 1 dO f dn de dm\ dfde__^ 
di ndi\ di'^dt 'dtj te dt'^ ' 



therefore -jz = — t-j~ + f 1 — i- j-jT 

dt dt \ h J dt 



nar 



(l^>(^-)S- 



C. p. T. 
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As in Art. 37, we may omit the term # -i-, if we bear in 
mind that the mean longitude will then be denoted by 
|wd^+€. Thus, on substituting the values of -r? and -^ , 
we obtain 

._ j.,„(«-„)h- 



■ -^(^;)(?-')<^-^}t- 



124. We shall now express the results of the preceding 
Articles in a form convenient for application ; and, for sim- 
plicity, terms of the order of the square of the disturbing 
force will be neglected. 

If u denote the excentric anomaly, we have 

r = a(l — 6Cosw) *.... (1), 

nt ■{- € " m =^ u — e ^mu (2). 



Hence (S)='^-^ 

aW — e COS u J 

fi 1 + g cos ^ 
a 1 — e cos w ' 

,, n ds /fl + ecosu\ . I , 

therefore j^ = "«/^( i_,co8 J ' "°'* "" "^''' 

And to transform the independent variable from t to u, 
we have from equation (2) 

dt 1 ,- . 1 dtc 1 dtsr 

-r^ =- (1-ecosw) T- + ~T7" y 

* du n ^ n au n au 
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of which the two last terms, being of the order of the dis- 
turbing force, may be omitted, since the terms in which the 
substitution is to be made are themselves of the first order. 

To obtain an expression for sin {0 — -or), which occurs in 
the formulae for the longitudes of perihelion and of the 
epoch, we have from equation (1), omitting terms of th$ 
order of the disturbing force, 



dr . du 



and from the equation 

- = ^ {1 + e cos (^ - «r)} (Art. 22), 

dv 
Equating the two values of -^ , 

. ,^ . hxi , du 
sm(& — -81)= — smtt-^. 

Hence, by substitution, the formulae of the preceding 
Articles become 

da 
du 



ia «7 8/:, V //l+«COStt\ 

^ = -2i/,a (l+ecos«)^^^-^— ^j. 
de -, ,- ,. //I + e cos mN 

d'GT _ 2^pa V(l — e') . / / I + e cos ^ \ 

dxjb" e \ \\ — e cos u) ' 

de 2kpa , ^ ,.^ ^. . . . ■ //I + e cos u\ 
du e ^ ^ ^ ^ V \1 — ecosW 

9—2 



•*^-^,-^_.^-w^_. ^.^--^ .^-^ " ' '- -- *^'*^ 
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125. These formubB are sufficient to determine the 
elemeuts of the orbit at any time, and being perfectly gene- 
ral, are applicable as well to the motion of comets, as to 
that of planets, but before we can integrate them, we shall 
require a knowledge of the form of p. Now the analogy of 
the terrestrial atmosphere would lead us to suppose that if 
the Sun be surrounded by an ethereal medium, its density 
decreases as the distance from the Sun increases. Moreover, 
the researches of Professor Encke on the comet which bears 
his name, seem to indicate the law of the inverse square. 
We will, however, merely assume p to be such a function 

of r, that when multiplied by a /(^j ) , and deve- 
loped in a series of cosines of t^ and its multiples, it takes the 
form 

j4 + 5^ cos ti + Cfe* cos 2ti+ ..• 
Thus our formulae become 
da 



du 

de 
du 



= - 2A;a' {A + liA+B)ecoau + ...}, 
= — 2ka\Acosu + -^ (l + cos2w) + ...l , 
«^=-2A;a-j^sinw + — sin 2u + .... I , 

-J- =ka(Aesia u, + ...). 

126. Supposing the orbit nearly circular, to examine tlie 
effects of the medium upon the elements of the orbit. 

Since the orbit is nearly circular, we shall neglect squares 
and higher powers of e ; thus the preceding formulae give 
on integration 



RESISTING MEDIUM, 133 

a = const. — ika^ {Au + {A+B)e sin w}, 

e = const. — 2A;a -|j.sint^ + -^ \u-\ ^ — jL 

enr = const. + 2ka -^-4 cos w+ -j- cos 2u\ , 

= const, — kaAe cos m, 

Hence in an entire revolution of the planet, the mean 
distance is diminished by 4i7rka*A, and the excentricity by 
iirkaBey while the longitudes of perihelion, and of the epoch, 

remain unchanged. Also from the formula n = ^, it ap- 

pears that the mean motion is, in an entire revolution, in- 
creased by GirknaA. 

127. We have already remarked that no traces of a 
resisting medium have yet been discovered in the motion of 
the planets : but, since k varies inversely as the ms^s of the 
body acted upon, the formulae of Art. 124 shew that such 
a medium, though too rare to influence the planets, might yet 
sensibly aflfect the motions of comets, in consequence of the 
extreme smallness of their masses. 



« •» 
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PKOBLEMS. 



1. Supposing in the Problem of the Three Bodies the 
relative orbit of two of the bodies to be a circle described 
uniformly, obtain equations for determining the motion of the 
third body; and transform the system of co-ordinates, so that 
the plane of the circular orbit being that of xy, the axis of x 
shall always pass through the two bodies in that plane. 

2. Shew that the plane of the orbit of a planet revolves 
about the planet's radius vector as an instantaneous axis*. 

3. A particle is describing an orbit round a centre of 
force which is any function of the distance, and is acted upon 
by a disturbing force which is always perpendicular to the 
plane of the instantaneous orbit, and inversely proportional 
to the distance of the body from the centre of the principal 
force. Prove that the plane of the instantaneous orbit re- 
volves uniformly round its instantaneous axis. 

4. A particle, acted on by a force varying as the inverse 
square of the distance and always tending to a fixed point, 
suffers slight disturbances : prove that generally there is a 
conic with the fixed point for focus, with which the body's 
actual path has a closer contact than with the instantaneous 
conic. 

* In this and the following problem, the plane of the orbit must be 
supposed to haye no angnlar veloeity about a normal to itself. See note to 
Art. 19. 
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5. Find whea the curvature of the instantaneous orbit 
of a body, acted on by disturbing forces, is the same as that 
of the actual orbit ; and shew that this is always the case 
when the only disturbing force arises from the action of a 
resisting medium. 

6. If R be expressed on the one hand as a function of 
r^, 5j, and z (Art. 11), and on the other as a function of r, 
0y i, and il, being measured on the plane of the orbit from 
the node, prove that 

dR_dR 

and obtain a formula for calculating the inclination. 

7. If i2 be expressed as a function of t and the usual 
elements, obtain the formulae 

dflt ^ 1 dR di_ coi(0^n) dR^ 
dt h sin i di ^ dt"^ h di ' 

where is measured on the plane of reference as far as the 
node, and thence on that of the orbit, and 

d0 ^ . otdn^ 



, .fdJ0 ^ . ^idll\ 



8. The central force being 4 + ^ , obtain the following 

equation for the apsidal motion, 

dl«r _ J a (1 ~ e*) ^' cos {0 — tsr) 
dt ehjfi r* 

a, e and 'sr being elements of the instantaneous ellipse. 

9. If the central force on a planet be -^ + ii\ shew that 
the planet may be supposed at any instant to be moving in an 
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ellipse, of which the meao distance is inversely proportional 

to the planet's distance from the Sun, provided ^=2pa, 

where p and a are two simultaneous values of r and a. Shew 
also that the latvs rectum of this ellipse is constant. 

10. When the disturbing function R is independent of 6, 
nnd expressions for -r. and -rr • 

If i2 = — , these expressions give variable values for e 

and -Bj, whereas the motion of the body actually takes place 
in a fixed ellipse: shew this, and explain the apparent 
paradox. 

11. A planet describes an orbit under the action of a 

force ^ tending to the Sun, /* not being quite constant: 

obtain the following equations for the variations of the ex- 
centricity and longitude of perihelion : 

If dfi be always positive, what in a whole revolution is 
the nature of its effects upon the excentricity and position of 
the major axis ? 

12. Prove that for any closed central orbit, the element 
introduced in the equation of vis viva is subject to periodic 
variations only, when the disturbing force is due to forces 
tending to centres and functions of the distances from those 
centres. 
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13. If the equation of the Moon's orbit be reduced to 
the form 

shew that the excentricity and longitude of perihelion may 
be found from the equations 

^ = -/sin(^-tsr), e^=/cos(^-w). 

Apply these equations to find e and -cr, when / is a small 
disturbing force, depending only upon the Moon's distance 
from the Earth, 

14. Assuming the differential equation for « in the Lunar 
Theory to be 

— + 5 = - -w«^ J- -I- - cos 2 {P'-mO)\ 

+ m^^jl sin 2 (^ - me^ , 

shew that if 7 be the longitude of the Moon's node, 

^ = — Tm* {1 — cos2 (m^ — 7) — cos 2 (^ — 7) 

+ co8 2(^-m^)}. 

From the above expression for ^, find the ratio of the 

mean motion of the node to that of the Moon, taking into 
account terms of the order «i*. 

15. If two planets disturbing one another were revolving 
in periods of 350 and 201 days, what form of terms in the 
disturbing function would demand examination ? 

16. The periods of Venus and the Earth are 224*7 and 
365 '256 days respectively; find approximately the period of 
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the long inequality anaing fix>ni their mutual perturbations, 
the important term in the disturbing function R being of the 
form 

PeV cos {13 {nt + e) - 8 {n't + e') - 3Br - 2^']. 

17. The radius vector of a planet is affected with a small 
periodical inequality; shew that its efifect may be represented 
by continued and periodical alterations of the excentricity and 

PT 

longitude of perihelion, the period of either being p — ^, 

where P is the period of the planet and T that of the in- 
equality. 

18. If in addition to the force of the Sun on a planet 
there be a small force tending towards the Sun, and varying 
inversely as the m^ power of the distance of the planet 
from the Sun, prove that the perihelion of the orbit will have 
a progressive or regressive motion according as m is greater 
or less than 2. 

Can you explain this result by reasoning similar to that 
used in Airy's Gravitation? 

19. It has been found by comparing theory with ob- 
servation that the perihelion of Mercury progresses at a rate 
greater by a than that due to the attraction of known bodies : 
shew that this increment would be accounted for if the law of 

force tending to the Sun were -a + ^> and if /a' = ac* a /-, the 

orbit being supposed to be nearly a circle, and the mean dis- 
tance to be c. 

20. The central force acting on a body being 



g+/*'^W. 
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shew to terms inclusive of fi and the square of the eccen- 
tricity, that the motion is in an ellipse revolving uniformly 
about the focus. 

21. Shew by means of the formula 

da _ 2na* dR 
dt fjL de * 

that the chief perturbation of the axis major of the Moon's 
orbit may be expressed by the equation 

where n and n are the mean motions of the Moon and Sun 
respectively. 

22. A satellite revolving in an ellipse of small excen- 
tricity is disturbed by another satellite revolving about the 
same primary; find approximately the variation of the mean 
distance and the motion of the apse, corresponding to the 
terms 

^V[l + 3cos {2 (n - 7iV+ 6 - e'}] 

in the function -B, having given 

da __ 2nd^ dR disr _ na \/(l — e*) dR 
dt /JL de * dt fie de ' 

23. Prove that, neglecting periodical variations, the ex- 
centricity of any orbit can always be represented by the 
diagonal of a parallelogram, whose sides are constant, and 
angle varies uniformly. ' 

24. Given the equations 

tan« { = N^' + jy/ + 2iV;ir, cos {h^t + S, - 8,), , 

^^""^ "" i\^;^\e + 5,)-hiV.cosS, ' 
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explain the nature of the motion of the node, when the mini- 
mum inclination is zero. 

25. Prove that as &.r as secular variations only are con- 
cerned the function F is constant, . . 

26. Considering only secular variations, obtain the fol- 
lowing equations : 

S ( m Va e* -ir) = C, Sf m hja tan* ^ "37 ) = C^« 

27. If the squares of the masses of two mutually dis- 
turbing planets were to each other inversely as their mean 
distances, shew that the nodes would oscillate through equal 
angles. 

' 28. If M, m, m' be the masses of three bodies mutually 
attracting according to the law of gravity, M being much 
larger than m or m\ and if v, v' be the velocities of m, m at 
distances r, r from the centre of if, supposed fixed, shew that 
the equation of vis viva for this case may be assumed to be 

2a and 2a being the major axes of the instantaneous ellipses 
of m and m\ 

29. Infer from the foregoing equation by the method of 
the variation of parameters the ratio of simultaneous changes 
in the mean distances and mean motions of two planets 
mutually disturbing. 

30. If r be the true radius vector, 6^ the projected lon- 
gitude, and X the latitude of a planet, obtain the following 
equation of motion r 

^-rcos«x(^j-ry +^=-^-. 
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31. Obtain the following equation between the pertur- 
bations of a planet in longitude and radius vector, whatever 
be the law of force, provided it be central and a function of 
the distance only, and provided such a function as R can be 
found : 

-4!FBr-2r~Br, 

ar 

where F denotes the central force, and h twice the sectorial 
area described by the undisturbed planet round the Sun. 

32. If the orbits of two planets which disturb each other 
be very nearly circular, shew that the inequalities of the 
radius vector may be immediately deduced from those of the 
longitude by means of the equation 

d.B0 

dt 

33. Integrate the equation 

^^^ + n\rBr==^^{Pcos{pnt+Q)}, 

determining the arbitrary constants so that Sr = 0, and 

d.Sr 

-^-— = 0, when t = : and shew that for small values of t, 

'Ff 



rSr = S (^ cos (3) , 



the case of i? = 1 being included. 
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34. A planet moves in a resisting medium of which the 
resistance 

J- (^\ . 

r'Kdt) ' 
apply the equation 

to obtain the following, in which e' is neglected : 
cP (rBr) ^ ^, ^g^ + n\r^.Se cos {nt + e- m) 



dR ^ 
dr 



df 



11 

+ 2/n*a {n< + -^ c sin (/if + € — «r)} = 0. 



35. The co-ordinates of the position at any time f of a 
disturbed planet being a? + S^, y+^y, ;^H- S^j, reckoned from 
the Sun's centre as a fixed origin, and referred to the plane 
of motion at a given epoch; and r being the heliocentric 
distance, x, y the co-ordinates of the position which the 
planet would have had at the time t, if the disturbance had 
ceased at the given epoch ; obtain the following equations for 
determining Sa?, Sy, Zz to the first order of the disturbing 
force ; 



df ' r' 



<? (Sz) , fiBz , dR' „ 



+ — = 
^ dy' "' 



in which /j, is the sum of the masses of the Sun and planet, 
and a^ is put for 

~m'ix'*-¥y'' + zTi + m''{Qc-tcy + (y-yy + (z-/y\-i, 

w! being the mass, and a?', y\ z' the heliocentric co-ordinates 
of the disturbing planet. 
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36. Shew that the eflfect of a resisting medium on the 
instantaneous orbit of a planet, would be to make the apsidal 
line regrede or progrede, according as the planet moved from 
perihelion to aphelion, or from aphelion to perihelion. 

37. Two small planets P, Q, very near each other, re- 
volve about the Sun in orbits very nearly circular, and make 
two revolutions about each other while they make one revo- 
lution about the Sun. Compare the sum of their masses with 
the mass of the Sun. 

If the line PQ move parallel to itself, what inference do 
you draw? 

38. If the motion of a planet round the Sun be disturbed 
by the action of another planet, the latter being supposed to 
describe a circular orbit of radius a with uniform velocity n, 
obtain the following exact equation : 



fdrX . , fd0X . fdzy « 2d0, , n 



2a 2rrbr 2m 

cos fi) + 



r a'* (a'« - 2aV cos w + r*)* ' 

where r is the radius vector of the disturbed planet, r^, 0^, z 
its co-ordinates referred to a fixed plane, and o) the inclina- 
tion of the radii vectores of the disturbed and disturbing 
planets to each other. 

39. Prove that, if the periodic times of a disturbed and 
disturbing planet are not commensurable, the secular changes 
of the orbit of the disturbed planet are the same as they 
would be if the mass of the disturbing planet were distributed 
over its orbit, in such a manner, that the part of the mass 
distributed over each portion of the orbit should be propor- 
tional to the time which the planet actually takes to describe 
that portion. 
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Additional Examples. 

40. If the disturbing force be radial and equal to fir^ 
prove that the mean distance at any time is 

a — 2a — ^ J^ cos (n< + € — -Bj) +7- cos 2 (n^ + € — -Bj)!-. 

{^Matli. Tripos, 1872.) 

41. A particle moves in a plane subject to a force directed 
to a fixed point in the plane, and depending on the distance 
from the point : prove these equations : 

fc'de hdb^^dR 
a» dt '^b^dt'^ da' 



h*de 
b'dt 


h da 

Vdt' 


dR 

- db' 


de 

"dt- 


-(i+e) 


dt "' 



where a and b are the intercepts of the direction of motion on 
fixed axes through the fixed point, A is an absolute constant, 
R is the potential of the force and e a quantity such that, t; 
being the velocity of the particle at time t^ v (< + e) is its 
distance measured along the direction of motion from the 
axis of a?. (M.T. 1873.) 

42. A large number of meteors are distributed uni- 
formly in the circumference of a circle of radius b concentric 
with the sun and in the plane of the planet's orbit ; shew 
that the perturbation in the mean distance of the planet is 
given by 

^ 2Ma /, o 6'\ / , . X 

oa = (l + f-jjj 6C0s(n^+6 — 'or), 

M being the whole mass of the meteors, a, e, n, e, -or being 
elements of the planet's orbit, and 6' : a', M* : fi^, and e^ being 
neglected. (M. T. 1880.) 
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43. If a, h be the semiaxes of the instantaneous ellipse 
of a planet's orbit, and the only disturbing force be a 
tangential one whose acceleration is^ prove that 

2a— — - 6 f— V = ^a*&. 
dt dt \dtj fjL 

(Jf. T. 1881.) 

44. If H be proportional to the excess of the moment 
of momentum of the orbital motion of a disturbed planet 
about the Sun in a circular orbit at distance a above the 
value of the same thing in the actual elliptic orbit ; then 
shew, from the equations for the variations of major-axis 

fJH r77? 

and excentricity, that, k being a constant, -j— = — A; -^— * 

at aw 

Also assuming that, when the inclination of the orbit to the 

plane of reference is zero, the variational equation for w is 

<?©• ^naj\— ^ dU 
dt fie de^ 



shew that -77- = k ttt * 

dt dU 



{M. T. 1882.) 



45. Explain how Elliptic Integrals and Bessel's Func- 
tions arise in the development of the disturbing function. 

(if. T. 1879.) 

46. If m, m denote the masses of two planets, 71, n\ 
their mean motions, and a, a', their mean distances from the 
Sun, and if 

m' —^G cos {{in ^ in) < + 7}, 

m—i(f cos {{i'n' ^ in) t + 7}, 
< a 

C. p. T. 10 
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represent any two corresponding terms in the respective 
developments of the portions 

m' -7j cos (r, r') and m -^ cos (r, /), 



of the two disturbing functions, prove that 

{SmiiKs Frizes, 1880.) 

47. Shew that the use of the ordinary disturbing func- 
tion is attended with certain inconveniences in the case 
where we have to determine the perturbations of a superior 
planet produced by the action of an inferior one. Also 
shew that these inconveniences may be avoided by referring 
the motion of the superior planet to the common centre of 
gravity of the Sun and inferior planet as origin, and by taking 

a quantity proportional to "w + - ~ as the disturbing 

function, where M is the Sun's mass, /a that of the disturbing 
planet, and jB, />, and r are the distances of the disturbed planet 
from the Sun, from the disturbing planet, and from the centre 
of gravity of these two bodies respectively. How would you 
develope the above disturbing function? (S. P. 1880.) 

48. Considering the secular variations of the excentrici- 
ties and perihelion longitudes of two mutually disturbing 
planets ; prove that in the case where m J a = mf Ja\ e and e' 
have the same greatest and least values E^ and E^, and «r -"bj' 

is never greater than tan"^ ( 9jp tp^ ' ) • 

(Jf. T. 1874.) 
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49. If there were only two planets moving about the 
Sun and at any instant the orbit of one were circular, the 
eccentricity of its orbit would always be given by 

c . cos (at + y3) 

where c, a, /3 are constants, and the longitude of its peri- 
helion would increase uniformly with the time. 

{M. T. 1876.) 

50. If 7 be the angle between the orbits of two mutually 
disturbing planets, prove that approximately, 



/ 9 . ' /— '2 . 4mm Jaa . sm' ^7 

is constant. (Pont^coulant, Bystlme du Monde^ Tome i. pp. 
458—460.) 

51. If I be the mean anomaly, v the true anomaly, g the 
angular distance of the perihelion from the ascending node, 
h the longitude of the ascending node in the fixed plane of 
reference, 

L — Jfjbay G = >Jfia{l'-' e^), H= Jfia (l—e!^) cost, 

the other letters having their usual meanings, and if li denote 

the ordinary disturbing function plus the term ^^, prove 

that a?, y, z, the co-ordinates of the planet referred to the fixed 
plane as the plane of osy and to the line from which the 
longitudes are measured as the axis of oc, are given by the 
equations 

Ox=: Grcos (y+g) cos A — ifr sin (v + ^r) sin A, 

Gy = 0r cos (v + g) sin h + Er sin {v + g) cos A, 

Oz « JO'-S'.rsm (v + g), 

10—2 
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and taking I, g, h, L, G, H, as the six elements, prove that 

dL^dB dG^dR dH_dR 
dt dl' dt dg' dt^ dh* 

{M. T. 1878.) 

52. If JIf, m be the masses of Sun and disturbed planet, 
€i the disturbing function due to the action of another planet; 
w, V the longitudes of the perihelion and ascending node, e 
the time of passage through perihelion ; a, e, i, the semiaxis 
major, excentricity and inclination of the orbit to a fixed 
plane ; further if 

" 2^' 

K = jM+m Ja (1 - e"), 

\==JM+m.Ja{l — e^) (1 —cost): 

then the pairs of elements (e, fi), {k, w), (v, X) will be 
canonical, i.e. 



de _ dCl dfi __ dfl « 
dt dfi * dt de ' 



(8. P. 1881.) 



I I »,l rij»i;^i j« Tltiijii jxi, 



APPENDIX. 

ON THE FOEM OF THE EQUATIONS OF ART. 39. 

1. On referring to Art 39, it will be seen that the 
formulae which have been obtained for calculating the elements 
of the orbit involve only partial dififerential coefficients of JB 
with respect to these elements, multiplied by functions of the 
elements which do not contain the time explicitly. As it is 
to this circumstance that these formulae mainly owe their 
advantage (since it renders them fit for use as soon as the 
partial differential coefficients have been calculated), it will 
be interesting to shew d priori that whatever system of 
elements be adopted, the formulae for their calculation can 
always be arranged in this way. 

2. If the motion of the planet be referred to three rect- 
angular axes originating in the centre of gravity of the Sun, 
we have the equations of motion (see Art. 9), 

d?^V"d^ ^ ^' 

^ fty_dR 

df^ t»~ dy *" ^ '' 



S^z fiz _ dR . . 

7^ "^Zji—T: W- 



df ' r* dz 
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Let a, h, c, d, e,f be the six elements introduced by in- 
tegrating these equations when jB = 0, and for -,- , -^ , -, 

write x\y, z\ then x\ y' and z' can be expressed as functions 
of t and the elements ; hence 

di \ dt) ^ dd dt ' dh dt 
where in {-7-) the elements are supposed constant. 
If in equation (1) we put R equal to 0, we have 



doif fdx\ dx da dx dh 



\dt) 






, ^ dx' /dx\ dR 

therefore _ _ ^_ j = _ ; 



\, - dR dx' da dx' db 

therefore ^ = 3a^ + d6^ + 



and similar equations hold for -j- and t- . 

Now since jB is a function of x, y and Zy 

dR^dR^ dRdy dRdz 
da dx da dy da dz da 

_ (dx do^ dy d'tf dz dz\ da 
\dai da da da da da) dt 

/dx dx' dy dy' dz dz\ ^ 
\da db da db da db) dt 

/dx dx' dy dy' dz dz\ dc 
\da dc da dc da dc) dt 

+ ;.. 
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3. "We may eliminate -^ from this expression: for, sup- 
posing X, y, and z expressed as functions of t and the elements, 
we have 

dx _ /dx\ dx da , ^^ ^^ , 
di'^Kdtj'^diL'dt "^db di "^ •'" ' 

but by the principles of the method of thie Variation of 

Parameters 

dx _ fdx\ 

di'^Kdi)' 

^, - dxda , dxdb dxdc , ^ 

therefore -r -ji + zii zj7 + J- zr*"^ "'^^' 

da dt do dt dc dt 

^. ., , dy da ^ dy dh dy dc ^ ^ 

dz da dz dh ,^^^1 _ n 
da dt db dt dc dt 

dx dii 
Multiplying these equations by -j- , ^- , &c., and adding, 

we obtain 

dx dx dy d]/^ dz dz\ da 



Ht 



+ 



da da da da da) dt 

fdx dx dy dy' dz^ dz\ db 

\db da db da db da) dt 

(dx dx' dy dy dz dz\ dc 

dc da dc da dc da) dt 



If this expression be subtracted from that for -7- in Art. 2, 
the latter may be written 

dR r rn ^^ . r n ^^ . 
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'"^L-'^J-dadft- 


dbda 


da db db da 

dz ds^ dz dz' 
da db db da ' 


Similarly, 






^^-rh 


„-\^j. 


rh /.1^^-L j&./> 



db "-'"^^dt ' L-» -J d« -•••' — 

4. By successive elimination between these equations, 
we can obtain expressions for -5-, -1-, &c., in terms of 

-Y-, -jj- y &c., [a, b\ [a, c], &c.: if, then, we can shew that 

[a, J], [a, c], &c., are independent of the time explicitly, it will 

follow that this is also the case with the coeflScients of -r- , 

dit D • xi- • ^ da db 

-yr, &c., m the expressions for -j- , -j-, &c. 
do at at 

5. To shew that [a, b] is independent of the time ex- 
plicitly. 

Let F= - ; then the equations of motion give 

(dx\^dy f^A^dZ. (^\=^ 
[dtj dx' \dt) dy* \dt) dz' ^ 

Now differentiating with respect to t only so far as it 
occurs explicitly. 



dp r -I _ <^ ^ (dx\ dx' d /dx\ 
dt^^' ^^"^ d^dt\db) '^'dbdtWJ 

^ dx dfdx\ dod d (dx\ 
db dt\da) da dt \db) 



+ 
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^dx d (dx\ Ax' dx' 
^ dadb\dtj db da 

^dx d /daf\ ^ ds^ ds^ 
dbda\dt/ da db 

+ 

^dxd^/dV\_dxd^/dV\ 
" dadb\dx) dbda\dx) 

,dyd^ /dV\ _ dyd^ (dV\ 
dadb\dy) dbda\dy) 

dz d (dV\ dz d fdV\ 
dadb\dz) dbda\dz) 

^dxd^ /dV\ dyd /dF\ dz d_ /dF\ 

"^ dadx\db) dady\dbj dadz\dbj 

_dxd_ /dV\ _d^d_ fdV\ _dzd^ (dV\ 

dbdx\ da) db dy\da) db dz \ da) 

da db db da 
Hence [a, J] does not contain the time explicitly. The 
same is of course true of [a, c], [J, c], &c. It follows, then, 
that whatever system of elements be adopted, we can always 
express their differential coefficients in terms of the partial 
differential coefficients of R with respect to them, multiplied 
by functions of the elements which do not involve the time 
explicitly, 

CoROLLABY. We now see that in the expression for [a, 5], 
viz. : — 

fdxdx' ^dxdx\ (dydy' dydy\ fdz dz ^ dz dz\ 
\da db dbda) \da db db da) \dadb dbda)* 
we may replace a?, y, z, x\ y\ z\ by x^, y^, z^, <, y/, <, their 
values at the origin of the time, without altering the value of 
this expression. We may next adopt instead of the constants 
a, 6, c, dy e, /, the same quantities x^, y^, z^, <, y^', z^. In 
this case we have evidently 
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K'> a-J = - 1, !>.', yj=-i, [<, ^o] = -i» 

and all the analogous expressions different from these will 
vanish. The final equation of Art. 3 of this Appendix will 
give 

dx^ _ <m dyo _ dR dz^ _ dR 

'dt^'dx;* It^dy^' Ht^dz^' 

^^_dR dyl^^dR dz^ ^ dR 
dt dx^ * dt "" dy^ * dt ~" dz^ * 

and no other combinations of the arbitrary constants will give 
simpler expressions than these for the variations of the 
constants adopted. From them the variations of the usual 
elliptic elements might be derived by proper analytical trans- 
formations. The group of constants aj„, y^, z^, x^, y^y z^ may 
be termed canonical elements. 

6. From the formula of Art. 3 of this Appendix, which 
is due to Lagrange, those of Chapter ii. may be deduced : for 
this we refer to Pont^coulant's Syst^me du Monde, Tome i. 
p. 542. 

Thus, if we denote by u the excentric anomaly of the 
planet, hy nt + I its mean longitude at the time t, by e the 
excentricity of the orbit, and by a the semi-major axis, con- 
nected with the constant n by the equation n^a* = fi, we shall 
have by the formulas of elliptic motion 

nt + l = u--esiiLU. 

Next, if we denote by X, F, the two rectangular co- 
ordinates of m referred to the plane of the orbit, the major- 
axis and latus-rectum being the co-ordinate axes, and the 
origin consequently at the focus of the curve, we shall have 

X = a cos u — ae, F= ajl — e^ sin u. 

Lastly, if we denote by v the angle which the major-axis 
of the orbit makes with the node or line of intersection of the 
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plane of this orbit and the fixed plane of xy, by H the angle 
which the node makes with the axis of a?, and by i the in- 
clination of the orbit to the plane of a^ ; we shall then have 
by transformation of co-ordinates, 

X = (X cos V — Fsin v) cos fl — (X sin i; + Fcos v) cos i . sin H, 
y = (Xcos 1/ — Ysin v) sin fl + (X sin v + Fcos v) cos i . cos fl, 
-ar = (X sin i; + Fcos v) sin i. 

These values of a?, y, 2? are very convenient for use in de- 
termining the values of the expressions [a, 6], [a, c], [J^ c], &c., 
since the constants in the expressions for aj, y, ^, are separated 
into two groups ; the variables X, Y only containing a, e, i, 
which depend on the form of the orbit and on the position of 
the planet on that curve at a given instant, being independent 
of the three constants v, H, %, which fix the position of the 
major-axis and of the plane of the orbit. 

By the aid of these values and of their partial diflferential 
coefficients substituted in forms such as 

fdxdx' dxdx\ fdydy' ^ dydj/X fdz dz_ dz dz\ 
\da db db da) \da db db da) \da db db da) 

we shall determine the values of the fifteen symbols [a, ZJ, 
[a, e\ &c., and by an easy calculation obtain the following 
values : 



r n aw r T J a (1 — e') ^ -, 



ae 



- ^_ COS % J a (1 — e'^) - ^- oecosi 

[t, n] = sin <f> Ja{l-- e% 

The values of the nine other symbols [a, e], [a, t], [e, Z], 
[^, fl], p, %\y \l, fi], [fl, J/], [Z, I/], [i, 1/] are equal to zero. 

After substituting these values in the general formula at 
the end of Art. 3 of this Appendix we shall find 
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dR __ andl J a (1 — ^) dp cos ija (1 — e') dfl 
da Ydi 2a di 2a dT' 

rfB_ a^ dv aecosi dn 

dR __ an da 
~dL^~2'dt' 

dR _ J a (1 — 6^^) da ae de 

dv 2a di " Ja (1 - e*) 5^ ' 

d-R cosi>ya(l — e^)da aecosi de . , /— 7= :2^ 

Whence, by an easy elimination, and observing that 
nV = /I, we deduce 

da _ 27ia' c?J2 
dt fji dl ^ 

de _an{l — e^) dR an J {I — e') dR 
dt e/jL dl eiM dv ' 

dl ^ 2a\ dR an ( 1 - e') dR 
dt jj, da fie de * 

dv _ anjl '— e* dJB an cos t d^ 

dt " efji de ^ sin i V(l - e^) d* ' 
dfl _^ aTi diZ 

^ "" /I sin i V(l - e') di ' 
rft an cos i dR an dR 

" dt fi sin i /(l-e*) dv ^ sin t V(l - e^) d^ ' 

formula3 from which we can readily obtain those of Art. 39 of 
this work when we take into consideration that 

Z = e — -or, i/ = 'cr — fl. 

For by Arts. 11 and 13, JB may be expressed as 

^ (a, e, i, r, — Q>, il), 
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or as either of <f> (a, e, i, w^ + €-tsr, tsr-H, n)= JS^, 

<f> (a, e, i, nt + I, z/, fi) = R^, 
and consequcDtly 

^,^_^ dR, dR^_dR^ 
dw dl "^ dv ' de dl ' 

dR^_dR^ dR, 



dSl dfl dv ' 

andso ^^dR, ^i = ^. + ^ 

dl de * dz/ dtiT de ' 

giZ, ^ dR, dR, dR^ 

da dil ■*" d^ "*■ d€ ' 
by aid of which the needful changes may be effected. 

NOTES TO ARTS. 31 AND 33. 

7. The process of Arts. 31 and 33 may perhaps present 
some diflSculty to the student, arising from the fact that the 
method of determining the position of the planet by its polar 
co-ordinates on the plane of the orbit, i.e. on a plane moving 
with the planet, is not applicable to the geometrical purpose 
of defining the position of a point on the curve of reference, 
a curve which, although it passes through the planet at some 
given instant, has nothing to do with the planet's subsequent 
motion. 

Of course the position of any point may be determined 
by its polar co-ordinates on a plane passing through it and 
some fiixed point or origin, the inclination of this plane to 
some fixed plane through the origin, and the longitude of 
its node. But inasmuch as an infinite number of planes 
can be drawn through two given points, some further con- 
dition is necessary to regulate the motion of that on which 
the polar co-ordinates are measured. If we are concerned 
with the motion of the planet, the condition is that this 
plane shall contain the direction of the planet's motion at 
the instant under consideration, in which case it becomes 



m^W 
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the plane of the orbit : if we are concerned with any other 
point, then some diflferent condition is requisite. In apply- 
ing this method to the curve of reference, we are, of course, 
at liberty to assign any condition we please, but it will be 
convenient that it be such that the co-ordinates of that point 
of the curve of reference which coincides with the planet 
shall be identical with those of the planet itself: this object 
is attained by the two conditions imposed in Arts. 31 and 33. 



ON THE EQUATIONS OF MOTION OF A DISTURBED PLANET. 

8. The principles of the conservation of areas and of the 
vis viva being applicable to the motion of the Sun and any 
number of mutually disturbing planets, four first integrals 
of the equations of motion may be found. Although these 
integrals are of little use in determining the motions of the 
several bodies, which cannot be completely found except by 
methods of approximation, we shall shew how they may be 
obtained in the case of the Sun and two planets. 

The differential equations of motion of the two planets 
relative to the Sun are 
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where /> = {(«;'- x)" + (y' - yf + {z' - z)*} \ 

From these we find 



m 



( d*z d'y\ , ,/ .d'z' , cPy'\ 



m 



and 

-=(M + m + m') mm \^p - ->, j (y^' - y'^). 
Hence 

mv + m'y f cPz , , cPz\ mz + m'z' ( d^y , cPy\ 
Mi-m-^m \ df dt* J M+m+m \ df df J 

Integrating, we have 

Tny + rny [ dz . ,dz\ mz + m'z' ( dy tdy\ 
Jf+wi + m V dt dtj -Jf+w + m'V dt dtj' 

with two similar integrals. 

These equations may be written thus 

+ »».'|(y'-y)|(/-;.)-(/-«)|(y'-y)} + 0,=0, 

with two similar equations. 

For a first approximation, in consequence of the magni- 
tude of the Sun's mass by comparison with that of any of the 
planets, we may neglect the term multiplied by mm', and 
we have 
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m 



/ dz dy\ , ,f,d2f ,dy\ 



m 



m 



J dx' 



dt 

d^ 

dt) 






Now 



/ dy dx\ 

y'dt-yn) 

is the projection of the area 

7(if+m)a(l-6^), 

described in the instantaneous elliptic orbit, on the plane of 
xy, to which that orbit is inclined at an angle i. Hence 

(a; f - y §) = J{M+m)a{l-^ (1 + W *)-* 

and we have 

m J{M + m)a{l- e') (1 + tan' t)"* 

+ m V(if+w')a'(l-e'") (1 + tan*i')"* = C3, 
or, retaining only terms aflfected by M, 

m^/a(l-c')(l + tanS•)-* 

+ m' J a' (1 - e *) (1 + tan' O"* = constant. 

If we multiply the diflferential equations (1) of motion of 
the planet m by the factors 

dx , dx' 

d^ if+w+m 



2m 



^_ 



dt 



2m. 



'^di'^'^'dt 



WT 



m-^-in 



~i 



dz ,dz 

J ^ "77+ m ^rr 

o a-2r -. dt dt 

dt M+m + m 
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and multiply also the differential equations (2) of motion of 
the planet m by similar factors, and add the six products, we 
shall find 



9-m (dx^dy^dz^^\ 
"^ U^ df'^ dt df^dtdf) 



« , fdz* d^x dv' d^v' d£ d^z\ 
\ at . air at ar at dv / 



don , dx 



^2"^* 



+ m \ dt dp) 



M+m 
dz ■ , d/ 

f'-di'^'^l^ f dl^z^ ,d^z'\ 

M + m + m \ dr dv 1 

^,-. fm dr m' dr\ 

dt\J{a,'-a!y+(y'-yf + (z'-z)y 
or integratiBg 

•»{(S)*^(S)"-(S)'}-1(S)'-(f)'^(S)} 

-2if 1 — + -7 — , == const* 

C.P.T. 11 
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au equation which may be written 



•i-mm 



\\dt dtJ ^\dt dti'^Kdt dtj) 



= const., 

an equation containing the principle of the via viva. 
Pont^coulant, Tome L, p. 210. 



EXAMPI^ OF NUMKRIOAL CALCULATION. 

9. By way of illustration, we will here shew How to 
calculate the inequality in the earth's radius vector due to 
the term • 

mV^ cos 2 {(w - wO « + € - €'] 

in By when Jupiter is the disturbing planet. 

Following the method of Chapter vlL, and taking the 
result of the first approximation, we havje (Art. 110) 

or, writing for P, its value, the coefficient of this term 
becomes . : j 

m' (dG, , 2nC, 

. . : . J , ■ « J^ : «_ 



(Sn-2n'){n-2n)\da a(n-n')j\ 

or since n^a^^fi^M+m, if M denote the Sun's mass, we 
may write, to the first order of the disturbing force, 

rnf^ n\ [ , (?(7^ 2na'(7, ] 

"^ if {3n - 2tt') (^ - 2.w') f da'^ n-n]* 
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a more convenient expression; for since n occurs in the form 
of a ratio, it is no longer necessary to suppose it expressed in 
•circular measure. 

dC 
The values of C^ and -^-^ must, of course, be found by 

the methods of Chapter ii. : by slightly changing our nota- 
tion, however, we shall be able to avail ourselves of the 
calculation of M. Bouvard, as given in the third volume of 
the M^canique Cilestey and the third volume of Ponte- 
coulant's TMorie Analytique du Systhne du Monde. 

I, 

Ifc2 = a'Cjj, a = -7 ,wehave 

a 

f * a!' da.^ a' da" a'^ da ' 

thus the coefficient becomes 

,< 

M{Sn-2n'){n-'2n) \ doL '^n-n') * 

or, taking the Sun's mass and the Earth's mean distance as 
the units of mass and of distance^ 

(3n-2n')(»-2n') f da "^ n^n j * 
We take from the Systlme du Monde, Tome iii;*, 

m' = .^^.g n = 1295977"-35, n' = 109256"-29, 

log,, a =1-2838993, log,, c, = 2-4501721, 

* In the second edition of the first yolnme, published in 1856, more 
recent yalaes of m\ n and n' are given, but our object being simply to illus- 
trate the method of calculation, the above values have been retained. 

11—2 
: • • :• 

^ A. ^ ^ 



> to. 



164 PLANBTABY THEORY. 

and with the aid of a table of logarithms the numerical value 
of the coefficient may easily be* found* Thus, in consequence 
of the term in R we have beeii considering, we find in Br the 
term 

- -0000092122 cos 2 {(n - w') « + e - e'}, 

ON THE METHODS OF CALCULATING THE MASSES OF 

THE PLANETS. 

10. There are in general two methods of determining the 
masses of the planets ; either by observations on a satellite, 
when the planet is accompanied by a satellite ; or by compar- 
ing the inequalities produced in their motion by their mutual 
action, as deduced from observation, with the same ine- 
qualities calculated from theory. The secular vacations are 
best adapted to give the most exact results ; but these are not 
yet known with sufficient accuracy to allow of this use. We 
are therefore obliged to recur to the periodic variations, and, 
by combining a vast number of observations, gather from 
them the most probable results*. 

11. When the planet is accompanied by a satellite the 
formula for calculating its mass may be obtained as follows : 

Let M, m, m' be the masses of the Sun, the planet, and 
the satellite : P, P' the periodic times of the planet about 
the Sun, and the satellite about the planet ; a, a' the mean 
distances of the planet from the Sun, and the satellite fiom 
the planet. Then we have 

P^ 27ra* p,^ iira'* 



therefore 



m +mf _ PV' 

* Pont^coolant* Systhne du Monde, Tome ni. p. 840. 
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or approximately, 

M "PV 

This equation gives the mass of the planet when that of 
its satellite is known. If the latter be neglected, the formula 
becomes 

12. In the case of the Earth, this method is not suf- 
ficiently exact, but the following may be employed. The 
attraction of the Earth on a body at its surface, in the 

parallel of which the square of the sine of the latitude is ^ t 

is very nearly the same as if the Earth were condensed into 
its centre. (See Pratt's Figure of the Earth, Art. 89.) Let 

then sin*Z = Q, 5r = the Earth's attraction on a body at its 

surface in latitude Z, h the mean radius of the Earth, E the 
mass of the Earth, M the ms|.ss of the Sun, P the length of 
the year, and a the mean radius of the Earth's orbit. Then 

E „ 27ra* 



6*' * ^JM' 



therefore 






©•■ 



where - = sine of Sun's parallax. 

For the methods of calculating the mass of the Moon 
we refer to Pontt?coulant's Systhme du Monde, Tome iv. 
p. 651. 
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ON THE CONSTRUCTION OF ASTRONOMICAL TABLES. 

13. In the present state of Astronomy, the values of 
the masses and elements, upon which these t-ables depend, 
may be considered as approximately known, so that small 
corrections only are necessary : we propose, in this note, to 
shew briefly how these corrections are eflfected. 

A series of observations are taken in R.A. and N.P.D., 
separated by considerable intervals of time: from each of 
these the geocentric longitude and latitude are obtained, and 
thence the heliocentric longitude and latitude. These co- 
ordinates are also calculated for each period of observation 
by the methods of the preceding Chapters, employing the 
existing values of the masses and elements. The differences 
between the observed and calculated values will be due to 
errors in the latter, for we assume that our observations can 
•be depended upon. From these errors in longitude and 
latitude, the corresponding corrections in the values of the 
masses and elements are obtained as follows : 

Let denote th6 longitude of the planet as calculated 
from theory, ^ + S^ its value as given by observation, m', m", 
&c., the existing values of the masses of the disturbing 
planets, Sm\ Zm'\ &c., S», &, &t, he the requisite corrections 
to the values of the masses and elements. Then, retaining 
only the first power of the excentricity and of the disturbing 
force, we have by the methods of the preceding Chapters, 

5 = w< + e + 26sin(n^+e-t!r) + 7w'P' + m"P" + , 

where m'P', 7n!'F\ &c., are the terms due to the perturba- 
tions of the planets m', m", &c. It follows that 

S^ = ^S/i+^S« + ^Scr + ^Se + P'Sw' + P"W+ , 

an de dv de 



